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CHAPTER  1 


Introduction 


1.1.  Literature  Review  and  Goals 

The  extension  of  classical  detection  theory  to  the  case  of  distributed  sensors  is 
discussed  in  [1];  in  particular,  the  problem  of  constructing  decentralized  Bayesian 
hypothesis  testing  rules  is  considered.  In  [2],  the  optimal  data  fusion  structure  is 
developed,  when  the  global  decision  is  obtained  by  weighting  local  decisions  accord¬ 
ing  to  the  reliability  of  detectors  and  comparing  to  a  threshold.  In  that  paper  optimal 
fusion  rules  are  derived  when  the  decision  rules  per  individual  detector  are  known. 
Those  rules  are  expressed  in  terms  of  the  probability  of  false  alarm  and  the  probabil¬ 
ity  of  miss.  The  systems  considered  in  [1]  and  [2]  have  a  fusion  center  which  always 
requires  all  the  sensors  to  transmit  their  local  decisions.  But  in  cenain  applications, 
such  continuous  communication  may  not  be  desirable;  such  is  the  case  in  environ¬ 
ments  with  adversaries.  In  [3],  one  of  data  fusion  methods  in  distributed  networks  is 
to  apply  the  Neyman-Pearson  approach  to  find  all  of  the  optimal  decision  rules  at 
each  site  (or  detector).  The  optimal  threshold  for  the  system  using  those  optimal 
decision  rules  found  at  each  site  is  not  stated.  In  [4],  the  problem  of  optimal  data 
fusion  in  the  sense  of  the  Neyman-Pearson  test  is  considered;  uncertainty  regions  at 
the  detectors  are  considered,  but  this  information  is  used  to  enhance  the  decision  at 
the  data  fusion  center.  A  region  where  a  definite  decision  cannot  be  made  is  called 
an  uncertainty  region.  There  are  no  communications  between  sensors  when  the 


1 


observation  falls  in  the  confident  region.  Papastavrou  and  Athans  [5]  evaluated  a 
two-sensor  network,  consisting  of  a  primary  sensor  and  a  consulting  sensor  using 
team  strategy  method,  with  performance  criterion  of  the  probability  of  error.  They 
also  provide  numerical  results  for  varying  quality  of  observations  at  different  sensors 
and  a  priori  probabilities.  The  relationship  between  the  position  of  threshold  in  the 
primary  sensor  and  the  system  probability  of  error  is  not  clearly  stated. 

Through  this  study,  we  applied  team  decision  strategies  to  three  different  sensor 
systems  and  an  analysis  of  each  system  was  performed.  The  three  different  systems 
are  a  two-sensor-system  (2SS),  a  three-sensor-system  (3SS),  and  a  two/three-sensor- 
system  (2/3SS).  The  main  goals  of  this  study  are  to  identify  the  level  of  risk  which 
prohibits  communication  between  sensors,  to  obtain  the  analytical  expression  of 
optimal  global  decision  rules  for  each  system  considered,  to  investigate  the  behavior 
of  decision  thresholds  and  system  performance  for  a  given  communication  cost  (or 
risk),  and  to  compare  the  performance  of  each  system  through  numerical  evaluations 
and  system  simulations. 

1.2.  Overview  of  Chapters 

In  Chapter  1,  general  concepts  are  discussed.  The  team  strategy  method  in  deci¬ 
sion  processes  and  the  communication  cost  involved  in  the  team  strategy  are 
described  in  this  chapter.  The.  binary  hypothesis  environment  and  assumptions  made 
in  deriving  the  expected  system  cost  are  also  stated.  A  couple  of  examples,  concern¬ 
ing  the  interpretation  of  the  communication  cost  constant  in  real  system,  are  also 
presented. 
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In  Chapter  2,  a  two-sensor-system  (2SS)  similar  to  the  system  studied  by  Papas- 
tavrou  and  Athans  [5]  is  considered.  The  model  consists  of  a  host  sensor  (HS)  and  a 
slave  sensor  (SS).  The  model  used  the  team  strategy  method  in  making  final  deci¬ 
sions,  depending  upon  the  local  decisions  made  by  the  host  sensor.  The  expected 
system  cost  is  expressed  in  general  probabilistic  terms.  This  expression  is  numeri¬ 
cally  evaluated  based  on  the  assumption  of  Gaussian  observation. 

Adding  an  additional  sensor  to  the  system,  a  three-sensor-system  (3SS)  is 
treated  in  Chapter  3.  In  this  system  both  slave  sensors  return  their  binary  information 
to  the  host  sensor  when  a  request  of  information  is  made  by  the  host  sensor.  The 
analytical  expression  as  well  as  the  numerical  evaluation  of  the  expected  system  cost 
are  also  performed. 

Chapter  4  contains  an  analysis  of  a  two/three-sensor-system  (2/3SS).  This  sys¬ 
tem  can  be  considered  as  a  combination  of  2SS  and  3SS  since  2/3SS  switches  from 
2SS  to  3SS,  or  vice  versa,  depending  upon  the  local  decision  of  the  host  sensor.  Most 
of  the  2/3SS  model  criteria  are  the  same  as  in  the  previous  chapters.  Numerical 
evaluation  is  also  done  for  this  system. 

The  results  from  the  numerical  evaluation  of  the  expected  cost  of  each  system 
are  presented  in  Chapter  5.  The  data  are  available  in  both  tables  and  plots.  The  plots 
are  attached  at  the  end  of  Chapter  2,  3,  and  4.  Each  system  is  compared  to  other  sys¬ 
tems  based  on  the  results.  In  the  numerical  evaluation  of  the  expected  system  cost,  C, 
FORTRAN  programs  are  written  and  these  are  attached  in  Appendix  A,  B,  and  C  in 
order  of  2SS,  3SS,  and  2/3SS,  respectively. 
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In  Chapter  6,  simulation  results  of  the  systems  analyzed  in  Chapter  2,  3,  and  4 
are  presented.  Tables  and  plots  are  used  to  show  the  data  from  the  simulation.  A 
FORTRAN  program  is  also  written  to  carry  out  the  simulation.  The  program  is 
attached  in  Appendix  E. 

Finally,  an  overall  summary  of  this  study  and  the  conclusion  are  written  in 
Chapter  7. 

U.  Environment 

A  binary  hypothesis  environment.  Ho  and  Hi ,  is  considered.  Hq  indicates  that 
there  is  no  target  present.  Hi  indicates  that  a  target  is  present. 

1.4.  Team  Strategies 

» 

The  sensor  communications  occur  only  when  the  host  sensor  declares  lack  of 
confidence  in  its  local  observation.  When  a  slave  sensor  transmits  only  a  binary  deci¬ 
sion  to  the  host  sensor,  some  information  received  at  the  slave  sensor  may  be  lost,  but 
the  risk  of  interception  by  adversaries  is  then  reduced.  Examples  of  the  costs  or  the 
risks  in  real  systems  are  given  in  section  1.6.  The  final  threshold  in  the  host  sensor  is 
evaluated  using  the  binary  decisions  transmitted  from  the  slave  sensors  and  a  prior 
probabilities  of  the  hypothesis.  The  final  threshold  (FT),  then,  is  compared  against 
the  observation  at  the  host  sensor  to  make  the  final  decision.  In  other  words,  the  final 
decision  at  the  host  sensor  is  declared  by  using  its  local  analog  data  and  the  binary 
decisions  transmitted  from  the  slave  sensors. 
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The  team  strategy  allows  collaboration  of  sensors  in  the  distributed  (multiple) 

sensor  system.  The  differences  between  distributed  sensor  systems  which  use  the 

team  strategies  and  those  that  do  not  are; 

(1)  The  host  sensor  in  team  strategies  has  an  overall  control  of  communication 
between  the  host  sensor  (HS)  and  the  rest  of  the  sensors,  namely  the  slave  sen¬ 
sors  (SSs). 

(2)  All  the  sensors  including  HS  and  SSs  are  capable  of  making  their  own  local 
decisions  utilizing  observations  from  their  local  environment. 

(3)  The  host  sensor  carries  muldple  thresholds  which  divide  the  decision  space  into 
either  three  regions  (2  thresholds)  or  five  regions  (4  thresholds).  It  uses  them  to 
recruit  inputs  from  other  sensors  accordingly,  which  means  that  the  communica¬ 
tion  schemes  between  the  host  sensor  and  the  slave  sensor  are  determined, 
depending  upon  the  decisions  of  the  HS. 

(4)  The  systems  do  not  have  a  central  data  fusion  center.  The  host  sensor  is  capable 
of  making  the  final  decision  either  based  on  its  local  observations  only,  or 
through  communication  with  the  slave  sensors. 

1.5.  Assumptions 

(1)  The  observations  received  at  different  sensors  are  mutually  independent  condi¬ 
tioned  on  each  hypothesis. 


(2)  All  the  sensors  used  in  the  model  are  considered  identical  in  performance. 

(3)  The  influence  of  the  number  of  observations  available  at  each  sensor  is  ignored. 

1.6.  Communication  Cost  Constant  (CCC) 

There  may  be  many  ways  to  interpret  an  application  of  the  communication  cost 
constant  (CCC)  in  real  systems.  An  example  would  be  the  communication  between 
sensors  in  a  hostile  environment,  where  interception  is  possible.  Then,  the  communi¬ 
cation  cost  constant  can  be  interpreted  as  the  probability  of  interception,  for  example. 

The  other  way  to  interpret  the  communication  cost  constant  is  that  a  limitation 
of  bandwidth,  a  duration  of  time  delay  in  communication,  quality  of  information 
obtained  by  communication,  etc. 

This  study  can  be  applied  with  a  modification  when  the  environment  is  non- 
hostile  and  the  cost  is  known.  In  this  case,  the  communication  cost  constant  can 
represent  a  physical  value,  such  as  a  dollar  cost,  etc.  For  example,  if  there  is  an  allo¬ 
cated  asset  or  capital  for  the  communication  between  each  party,  the  asset  (or  budget) 
should  be  wisely  used  to  obtain  the  information  from  the  other  party.  If  the  asset  is 
SlOO.OO/month  and  the  communication  cost  is  $  10.00/communication,  there  are  only 
10  communications  per  month  allowed.  Thus  the  system  should  use  the  communica¬ 
tion  capability  when  it  is  really  required.  On  the  other  hand,  if  the  communication 
cost  constant  is  SI. 00/communication,  the  system  have  100  communications  per 
month.  This  case  the  system  can  use  the  communication  capability  more  frequently. 
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In  this  paper,  the  communication  cost  constant  is  interpreted  as  the  probability 
of  interception.  The  larger  communication  cost  constant  indicates  greater  risk  in 
communication  with  other  panies.  Thus  when  the  communication  cost  constant  is 
null,  the  communication  between  parties,  i.e.,  the  host  sensor  and  the  slave  sensors, 
are  encouraged  and  desirable;  however,  as  the  communication  cost  constant 
increases,  the  exchange  of  information  is  restricted  to  the  cases  of  "must  communi¬ 
cate"  only. 
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CHAPTER  2 


Analysis  of  a  T wo-Sensor-System  (2SS) 


2.1.  The  Model  and  Configuration 

2.1.1.  Host  Sensor  (HS)  and  its  Decision  Maker  (HDM) 

The  host  sensor  has  two  decision  boundaries  providing  three  decision  regions. 
One  of  the  boundaries  is  called  a  lower  threshold  (TL)  and  the  other  is  called  an 
upper  threshold  (TU).  When  an  observation  falls  below  TL,  HDM  will  declare  "No 
Target  Detected”.  When  an  observation  is  between  TL  and  TU,  HDM  declares  "Not 
sure  and  communication  necessary".  Finally,  when  an  observation  falls  above  TU, 
HD.M  declares  "Target  Detected”.  Throughout  this  paper,  the  three  decisions  men¬ 
tioned  in  the  above  will  be  denoted  by  a  set  Uhs  =  {0,  ?,  1},  respectively.  These 
decision  regions  are  shown  in  Figure  2.1.  These  thresholds  can  be  varied  from  one 
mission  to  another,  depending  on  the  specific  requirements  and  constraints.  The 
thresholds  control  decision  accuracies  and  the  frequency  of  communication  between 
the  host  sensor  and  the  slave  sensors.  The  narrower  the  gap  between  TL,  and  TU  is. 
the  less  communication  between  HS  and  SS  would  occur.  This  is  because  the  gap 
between  the  thresholds  is  directly  related  to  the  uncenainty  decision  region  in  the 
host  sensor.  The  decision  region  between  TL  and  TU  may  be  called  a  dubious  region 
or  uncenaintv  region. 
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TL2 


TU2 


Figure  2.1  Decision  Boundaries  of  HS  for  2SS 


2.1.:.  Slave  Sensor  (SS)  and  its  Decision  Maker  (SD.M) 

The  slave  sensors  have  a  single  decision  threshold  (Tss)  providing  two  decision 
regions.  The  slave  sensor  does  not  have  an  uncertainty  region,  meaning  the  SD.Ms 
are  forced  to  make  a  decision  either  "O  "  or  "1". 

WTen  an  observation  falls  below  Tss.  SDM  declares  "No  Target  Detected”. 
WTen  an  observation  falls  above  Tss.  SDM  declares  "Target  Detected".  In  this  paper 
these  decisions  are  represented  by  a  set  L'ss  =  (0.  1 }. 

2.1.3.  Host  Sensor's  Final  Decision  Threshold 

When  the  communication  between  the  host  sensor  and  the  slave  sensors  occurs, 
the  analog  data  of  the  host  sensor  and  the  slave  sensor’s  binary  data  are  used  to  deter¬ 
mine  the  final  decision.  The  threshold  for  the  f  nal  decision  is  evaluated  utilizing  the 
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binary  information  from  the  slave  sensors  and  a  priori  probabilities.  Then  this  final 
threshold  is  compared  against  the  observation  received  by  the  HS  to  determine 
whether  there  is  a  target  or  not.  The  final  decision  is  denoted  by  Up-  FT  can  vary 
from  one  evaluation  to  another  since  Uss  provided  from  the  SDMs  may  differ  from 
one  communication  to  the  other. 

2  1.4.  The  Overall  Process  of  The  Model 

An  observation,  yns.  which  is  received  at  the  host  sensor,  is  mutually  indepen¬ 
dent  from  the  observations  received  by  other  sensors.  When  the  observation  is 
greater  than  or  equal  to  TU  (TU2  in  Figure  2.1)  or  is  less  than  or  equal  to  TL  (TL2  in 
Figure  2.1),  the  host  sensor’s  local  decision,  Uhs.  i  or  0,  respectively,  becomes  the 
final  decision.  Up.  When  yns  is  between  TL2  and  TU2,  a  request  of  assistance  from 
the  host  sensor  to  the  slave  sensor  is  transmitted.  The  slave  sensor  returns  the  local 
decision,  Uss,  ^o  the  host  sensor’s  team  processing  unit  upon  the  request.  This  com¬ 
munication  process  involves  a  communication  cost  constant  (CCC).  L’ss  ts  also 
determined  based  upon  an  independent  observation  at  the  slave  sensor.  The  slave  sen¬ 
sor  makes  a  binary  decision  since  it  only  has  one  decision  threshold.  Refer  to  Figure 
2.2. 

2.2.  Definition  of  the  System  Cost  Function 

The  system  cost  is  defined  by  the  total  system  probability  of  error.  The  follow¬ 
ing  is  the  cost  function  of  the  system,  C(,.„),  which  is  represented  by  error  probabili¬ 
ties  of  the  individual  sensor  as  well  as  the  error  caused  by  the  team  process. 
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a 


Figure  2.2  Model  Configuration  of  2SS 


C(z,yHs.TL,TU,CHS:Ss)  =  d  -z)  •  Pens  +Cus:ss)  (2.2. 1) 


In  (2.2.1),  z  is  a  value  that  determines  whether  the  communication  should  be 
made  or  not.  z  takes  a  binary  number  either  0  or  1.  When  the  host  sensor  makes 
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decision,  z  becomes  1.  In  case  the  host  sensor  makes  a  confident  decision,  z  becomes 
0.  It  is  obvious  from  the  expected  system  cost  function  that  the  team  decision  opera¬ 
tion  takes  a  role  only  when  a  communication  channel  is  open,  i.e.  z  =  1,  between  the 
host  sensor  and  the  slave  sensor.  When  z  =  0,  meaning  that  there  will  be  no  commun¬ 
ication  between  the  host  sensor  and  the  slave  sensor,  the  cost  function  becomes  that 
of  a  centralized  system  of  single  sensor  with  a  possibility  of  smaller  error. 

23.  Evaluation  of  an  System’s  Expected  Cost,  C 

Since  the  system  cost  function  is  defined,  it  is  possible  to  evaluate  an  system’s 
expected  cost. 

C=E{C(z,TL,TU,Chs;ss)} 

=  C(0.TL,TU.Chs:Ss)  •  Pr(z=0)  +  C(I.TL,TU,CHS:SS)-Pr(z=l) 

=  C(0.TL,TU.Chs:Ss)-  {1-Pr(z=l))  +C(l.TL,TU.cc)  P,(z=l) 

=  C(0.TL.TU.Chs:Ss)+  { C(  1  .TL.TU,Chs:ss) -C(0,TL.TU.cc) )  ■  P,(z=l )  (2.3.1) 

By  evaluating  C(0„„)  and  C(l„„),  and  using  (2.2.1);  the  following  is  obtained; 

C(O.TL.TU.CHS:Ss)=PeHs  (2.3.2) 

C(l,TL,TU,CHs.ss)  =  PeT«a„,  +Chs:SS  (2.3.3) 

Therefore  the  expression  of  equation  becomes  as  follows: 

^  ~  Pens  +CHS;SS-PeHs) '  Pr(z=l)  (2.3.4) 

This  equation  is  further  developed  in  detail  such  as; 

Pens  =Pr(faJse  local  decision  at  HS) 
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=  Pr(DecideHi  I  He)  ■  Pr(Ho)  +  P,(Decide  Ho  I  Hi)  Pr(Hi) 

=Pr(yHS^TU  I  Ho)  •  P,(Ho)  +  P,(yHS^  I  Hi)  •  P,(H,)  (2.3.5) 

Pj(z=l)=p,(uncertainty  in  decision;  communication  channel  open) 

=Pr(TL<yHs<TU) 

=  P,(TL  <  yHs  <TU  I  Ho)  ■  P^Hq)  +  P,(TL  <  yws  <  TU  I  Hj  )  •  P,(Hi  )  (2.3.6) 

^^Teun  (error  resulted  by  communication  using  team  strategy)  =  Pr(E) 

=P,(E  I  yHs  e  [TL.TU1)  •  P,(yHs  e  [TL.TU]) 

=P,(E.  yHse  [TL.TU]  I  Hq) •  P,(Ho)  +  P,(E.  yHs  e  [TL.TU]  I  H,)  P,(Hi) 

=  P,(E.  yHse  [TL.TU]  I  Hq.  Uss)-Pr(Uss  '  Ho)-Pr(Ho) 

+  P,(E.  yHs  €  [TL.TU]  I  H, .  Ujs)  •  Pr(Uss  I  H, )  •  Pr(Hi )  (2.3.7) 

where  Pr(Ho)  and  Pr(Hi)  are  a  priori  probabilities  of  the  environment  Hq  and  Hi 
respectively.  Uss  is  the  local  decision  determined  by  the  slave  sensor. 

In  the  (2.3.7)  it  is  quite  reasonable  that  a  decision  of  the  slave  senior,  the  binary 
data,  takes  a  pan  since  the  communication  between  the  host  sensor  and  the  slave  sen¬ 
sor  is  established  as  a  team  effort.  This  is  shown  in  the  equation  by  giving  the  proba¬ 
bility  terms  conditioned  on  the  decision  of  the  slave  sensor.  To  evaluate 

P,(E,  yHs  e  [TL.TU]  I  Uss.H,).  i=0,l  (2.3.8) 

it  is  necessary  to  compute  the  likelihood  ratio.  A(yHs,Uss).  This  probability  is  the 

probability  of  error  induced  in  the  host  sensor  due  to  the  communication  with  the 

slave  sensor  (refer  to  (2.3.7)).  Thus,  the  probability  of  error  in  the  slave  sensor  will 

contribute  to  the  probability  of  error  evaluated  in  the  host  sensor. 
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2.4.  The  Likelihood  Ratio  Test 

In  evaluating  the  Likelihood  Ratio  (LR),  it  is  assumed  that  the  individual  obser¬ 
vations  received  at  each  sensors  are  independent  form  each  other.  Thus,  their  local 
decisions  are  also  independent.  In  other  words,  the  local  decision  of  the  slave  sensor 
is  statistically  independent  from  (or  not  coupled  with)  either  the  local  decision  or  the 
observation  of  the  host  sensor.  Then  the  LR  of  this  system  can  be  written  as 


A(yHs.Uss) 


PrCyHs.  Uss  I  Hi) 

Pr(yHS*  Uss  I  Ho) 


(2.4.1) 


and  the  above  equation  is  re-written  as  follows: 


A(yHs.Uss) 


Pr(yHS  I  Hi) 
Pr(yHS  1  Ho) 


Pr(UsS  I  Hi)  Pr(Hi) 
PrOJsS  I  Ho)  '  Pr(Ho) 


Up=l 

> 

< 

Uf=0 


where. 


Cio  ~Uoo 

Coi  -Cn 


a  pre-calculated  threshold,  and 


Cao  :  a  cost  of  deciding  a  given  (3 


(2.4.2) 


(2.4.3) 


For  the  most  of  cases  it  is  assumed  that  a  cost  of  making  a  false  decision  and 
cost  of  missing  target  are  the  same,  i.e. 

Qi  =Cio, 

and  this  also  applies  to  a  cost  of  making  a  correct  decision,  for  example, 

Coo  =  Cii. 

These  conditions  will  give  the  pre-calculated  threshold  Xi  =  l.  Re-writing  (2.4.2)  into 
(2.4.4),  which  plays  an  important  role  in  evaluating  (2.3.8)  is  obtained. 
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PrCyHS  I  Hi) 
PrCyHS  •  Ho) 


Pr(HQ)  Pr(UsS  I  Ho) 

P,(Hi)  ■  P,(Uss  I  Hi) 


(2.4.4) 


and  now  defining  the  following, 


g(yHs) = 


Pr(yHS  I  Hi) 
PrCyHS  •  Ho)  ’ 


f(Uss)  =  >.fTo 


PrCUsS  I  Ho) 
PrCUss  I  Hi) 


(2.4.5) 


(2.4.6) 


(2.4.7) 


then  (2.4.4)  becomes  as 


Uf=1 

g(yHs)  <  f(Uss) 
Uf=0 


(2.4.8) 


g(yHs)  is  the  decision-statistic  and  f(Uss)  is  depends  on  the  slave  sensor’s  decision. 
The  function  f(Uss)  represents  the  final  threshold  (FT)  in  the  host  sensor  after  a  com¬ 
munication  is  exchanged.  As  is  seen  in  (2.4.8),  the  function  f(Uss)  t^kes  two  dif¬ 
ferent  values  (thresholds)  depending  on  the  decision  of  the  slave  sensor,  Uss-  Since 
the  slave  sensor  is  forced  to  make  a  decision  based  only  upon  its  observation,  Uss  is 
going  to  be  either  "0”  or  "1".  More  explicit  expression  of  the  function  f(Uss)  is 


f(Uss=0)  =  kiTo 


f(Uss  =  l)  =  XfTo- 


P,(Uss=OIHo) 

P,(Uss=OIHi)’ 

P,(Uss=llHo) 

P,(Uss=llHi) 


(2.4.9) 


(2.4.10) 


Then,  the  final  decision,  Up,  rules  can  be  written  as 
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Uf=1 


1  .  if  Yhs  ^  TU 
1  .  if  g(yHs)  ^  f(Uss) 

0  .  if  g(yHs)  <  f(Uss) 

0  .  if  yHs  <  TL 

All  equations  needed  to  evaluate  (2.3.8)  which  is  from  (2.3.7)  are  obtained.  By 
evaluating  (2.3.8)  further,  the  following  is  derived: 

P,(E,  yHi  ^  [TL.TU]) 

Uss“t 

=  Pr(Ho)  •  s  Pr(E.  yns  e  [TL,  TU]  I  Uss .  Hq)  •  P,(Uss  1  Ho) 

Uss-o 

Ujs-l 

+  P,(H,)-  I  P,(E.yHs€[TL,TU]IUss.Hi)-P,(UsslH,) 

Uss*0 

uf=i 

=  P,(Ho)  •  X  Pr{g(yHs)  f(Uss)  and  yws  e  [TL,  TU]  I  Uss-Hq]  •  Pr(Uss  I  Hq) 

t-'ss=0 

L'sjal 

+  P,(H,)-  X  Pr{g(yHs)u‘Lof^ss)andyHs6  [TL.TU]  IUss.H,)-P,(Uss  IH,)(2.4.11) 
Thus,  when  (2.3.5),  (2.3.6),  (2.3.7),  and  (2.4.11)  are  substitute  into  (2.3.4),  the  gen¬ 
eral  expression  of  the  expected  system  cost  is  obtained. 

2.5.  C  under  Gaussian  Assumption 

The  Gaussian  distribution  are  applied  to  the  probabilistic  expression  of  C  so  that 
numerical  method  can  be  used  to  evaluate  the  C  for  various  thresholds  in  the  host 
sensor  that  effects  system  performances. 

2.5.1.  The  Q(y)-function 

In  evaluating  the  probabilities  which  are  involved  in  the  analytical  expression  of 
the  expected  system  cost,  an  integration  of  Gaussian  probability  density  function  is 
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required.  We  define  the  Q(y)  function  to  be 

am  *2^ 

Q(y)  =  J^e  2  dz  (2.5.1. 1) 

y  ^2)C 

2.5.2.  Probability  Density  Functions 

The  probability  density  functions  under  Hq  and  Hi  at  each  sensors  are  written 
below.  For  the  symbols  used  in  the  expression,  please  refer  to  the  beginning  of  the 
thesis  under  "Symbols  used  in  Chapter  2". 

2.S.2.I.  Gaussian  PDFs  at  the  Host  Sensor 


-(yHs-^HSpr 

fHSo(yHs)  =  -7=^^ - e 


HSo 


-(yHs-MHSi  r 


fHS,(yHs)  =  ^==^^ - e 


(2.5.2.1.1) 


(2.5.2.1.2) 


2.5.2.2.  Gaussian  PDFs  at  the  Slave  Sensor 


-<yss-Hsso) 

hso(yss)  =  -^ - e 


"'2Jcasso 


-Kyss-Mssi)^ 


fssi  (yss)= ~r— C 
"'2rtOsSo 


(2.5.2.2.1) 


(2.5.2, 2.2) 


2.5.3.  Decision  Boundary  of  Slave  Sensor  (SS) 

A  LRT  is  used  to  find  SD.VI’s  optimal  decision  threshold,  Tss- 
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,  ,  Pr(Ho)  Qo-Coo  Pr(Ho) 

ss(yss)  p,(Hi)  ■  Coi -Cn  '  Pr(Hi) 


(2.5.3. 1) 


and  the  LR  can  also  represented  as 


^ss(yss)  = 


Pr(yss  I  Hi) 
Pr(yss  I  Ho) 


1 

r  " 

(yss-itss,)^ 

f  cxp^ 

V2jcosSi 

2assi 

J 

(yss-)tsso)^ 

cxp^ 

V2jtosso 

2ctsso 

w  ^ 

(2.5.3.2) 


Equating  (2.5.3. 1)  and  (2.5. 3.2)  and  taking  the  natural  logarithm  on  both  sides  of  the 
equation,  we  obtain 


lOgej^Kss 


P,(Ho)] 

■  P  (Hi)J 


=  10ge 


*JSSo 

<JSSi 


(q|si  -q|so)yss^-*-2(ussi-cr§so-Psso-q|si)yss-*-n§so-g|si  -nlsi  -giso 

2oisooisi 


By  re-arranging  the  terms,  the  above  equation  can  be  written  as  (2. 5. 3. 3). 


(o5s  1  -  o§So  )yss  ‘ + 2(pss  I  '*3550  “  l^sso  '^^ssi  )yss 


Pr(Ho) 


H?So-<J^Si -^hi  -^Iso  -2aiso-aisiloge|Xt5j  •  =  0  (2.5.3.3) 

Solving  (2.5. 3. 3)  for  yss,  the  optimal  threshold  for  the  slave  sensor,  Tss  is  found. 


^ss  - 


■<Mssi  'oiso  -Psso'<j|si ) 
(o|si  -also) 


(a|s,-<j|so) 


V 


4iss,  olso  -Hsso  Ois,)*-(ols,  -olso) 


tiJso'ois,  -Risi-asso  -2o|so  als,  lod 


P,(Ho)T 


(2.5.3.4) 


In  (2.5 .3.4)  we  require  OsSq  equal  •  Por  the  case  Ossg  =  =  ct,  (2. 5. 3. 3) 


becomes 


fp  j 

\s  •  =2a^(^issi  -WSo)yss -M^So) 

Again,  solving  for  yss . 


l^SSo+l^-SSi 
Tss  = - r - +  • 


■iogci  ^tss  ■ 


Pr(Ho) 


(2.5.3.5) 


2  USSi-USSo  I  Pr(Hi)J 

For  the  numerical  evaluation  performed  later  in  this  chapter,  the  threshold  for  the 

slave  sensor  is  evaluated  by  using  (2.5. 3.5).  The  decision  at  the  slave  sensor  is  carried 

out  as  below: 


fo  •  if  yss  <Tss 

.ifyss^Ts,  C.5,3,6) 

2.5.4.  Calculation  of  C  for  Gaussian  Model 

Let’s  represent  the  equations  derived  in  the  previous  sections  using  Gaussian- 
distributed  data. 


PcHs  =Pr(yHs2:TU  I  Ho)-P,(Ho)  +  Pr(yHS^TLI  Hi)  P,(Hi) 

oo  TL 

=  Pr(Ho)-  J  fHSo(yHs)  tiyns  +  Pr(Hi)-  JfHSi  (yHs)tiyHS 
TU  — 


=  Pr(Ho)-Q 


tu-phso 


L 


<^HSo 


+  P,(H,)H  1-Q 


TL-tiHSi 


OHSi 


(2.5.4. 1) 


P,(z=  I )  =  P,(TL  <  Xhs  <  TU  I  Ho)  •  Pr(Ho)  +  Pr(TL  <  y„s  < TU  I  H, )  •  P,(H, ) 
TU  TU 

=  Pr(Ho)  •  J  fHSo(yHs)  tiyiis  +  Pr(Hl)  •  f  ThSi  (yus)  tiVHS 

TL  TL 
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=  Pr(Ho)- 


<• 

TL-UhSo 

TU-)iHSo 

V 

Q 

<3HSo 

OHSo 

+  Pr(Hl)i 


TL-^lHSi 

-Q 

TU-M-hs, 

V 

1  ^ 

Ohs, 

^^HSi 

J 

k 

(2.5.4.2) 


and 


‘  ®Tci 


U»»l 

=  P,(Ho)-  S  Pr{g(yHs)>-^'Uss)andyHse  UL.TU]  I  Uss.Hq}  -PtCUss  I  Ho) 

u„^ 

+  P,(Hi)-  2  P,{g(yHs)<f(Uss)andyHs6  [TL.TUl  I  Uss.Hi)-P,(Uss  I  Hi) 
Un»0 

«  Tss 

=  P,(,Ho)-  J  fHSo(yHs)dyHS‘  J  fss0<yss)dyss 

Wss-O)  — 


+  Pr(Ho)-  I  fHSo(yHS)dyHS'  J  fsSo(yss)dyss 

f(Uss  =  1)  Tss 

f(Uss»0)  Tss 

+  Pr(Hi)-  /  fHSi(yHs)<iyHs- J  fssi(yss)dyss 

f(Uss  -  1)  O* 

+  Pr(Hi)-  I  fHSi(yHs)dyHS  -  J  fssi(yss)dyss 

—  Tss 


f 


=  Pr(Ho)-Q 


f(Uss=0)-)lHSo 


<JHSo 


Q 

1- 

'1 

Tss  “M-sso 

osso 

J 

. 

k  > 

J 

•P,(Ho)Q 


f(Uss  =  l)-dHSo 

.  n 

Tss-M-sSo 

OHSo 

k 

V 

Osso 

k  « 

20 


f 

*  ' 

« 

f 

+  Pr(Hl)- 

1-Q 

f(Uss=0)-PHSi 

1-Q 

Tss-M-ssi 

y 

<7HS, 

<JSSi 

. 

b  • 

+  Pr(H,)- 


1-Q 


* 

•  « 

f(Uss  =  l)-dHSi 

'Q 

Tss~M-SSi 

<7HSi 

OSSi 

* 

(2.5.4.3) 


Substituting  (2.5.4. 1),  (2.5.4.2),  and  (2.4.5.3)  into  (2.3.4),  the  expected  system  cost  is 
obtained. 


2.6.  Numerical  Evaluation  of  C 

In  the  previous  section  C  is  represented  in  terms  of  Q(y)-function,  which  makes 
possible  to  evaluate  C  numerically.  The  purpose  of  the  numerical  evaluation  is  to 
determine  the  expected  system  cost  at  the  various  thresholds  position  in  the  host  sen¬ 
sor,  meaning  the  position  of  TL  and  TU,  so  that  the  optimal  thresholds  for  the  system 
can  be  realized  with  different  communication  cost  constants.  At  the  optimal  thres¬ 
holds  the  expected  system  cost  is  minimum. 

The  a  priori  probabilities  of  the  environment  are  considered  equiprobable, 
Pr(Ho)  =  Pr(Hi)  =  0.5.  For  the  statistics  of  the  observations  we  take 
*^HSo  =  ^HSi  =  *^SSo  =^SSi  =<7=1.  The  mean  values  of  the  observations  at  each 
sensor  are  |iHSo  =  l^sso  l^HSi  =4SSi  =  1-  The  communication  cost  con¬ 

stant  is  held  a  constant  value  until  all  the  expected  system  costs  are  evaluated  at  the 
desired  threshold  positions.  The  thresholds,  TL  and  TU.  are  varied  with  a  relation¬ 
ship  of  TU  =  -TL.  This  threshold  relationship  is  selected  because  of  the  symmetric 
nature  of  Gaussian  PDF  and  its  observations.  The  results  from  the  numerical 
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evaluation  of  C  are  plotted  in  Figure  2.3,  Figure  2.4,  Figure  2.5,  and  Figure  2.6.  For 
the  tabulated  data  of  Figure  2.5,  please  refer  to  Table  5.1  in  Chapter  5.  The  computer 
program  is  written  for  (2.3.4),  substituted  with  Pe^s*  Pr(z=l).  and  Pex„n,-  which  are 
expressed  in  (2.5.4. 1),  (2. 5.4.2),  and  (2.5. 4.3),  respectively.  The  program  is  listed  in 
Appendix  A. 

2.6.1.  Comments  on  Numerical  Evaluation 

The  system  expected  costs  are  evaluated  over  the  threshold  positions  on  the  host 
sensor’s  observation  space  for  a  given  communication  cost  constant.  The  following 
figures  are  plotted  from  the  results  obtained  through  the  numerical  evaluation  of  C. 

Figure  2.3  shows  the  expected  system  cost  as  the  threshold  is  departing  from  the 
origin  (position  0.0)  for  each  communication  cost  constant,  CCCl.  TU  moves  in  the 
positive  direction  and  TL  moves  in  the  negative  direction.  When  CCCl  >  0.5,  in  the 
curve  of  the  expected  system  cost  vs.  HS  threshold  position,  the  expected  system  cost 
never  gets  smaller  than  the  cost  at  threshold  position  0.0.  The  dotted  curve  indicates 
the  maximum  CCCl  which  has  a  minimum  other  than  at  the  threshold  position  of  0.0. 
Figure  2.4  is  the  enlarged  version  of  Figure  2.3,  where  the  minima  of  the  curves  are 
shown  clearly. 

Figure  2.5  is  a  plot  of  extracted  information  from  Figure  2.3.  It  shows  the 
behavior  of  the  minimum  expected  system  cost  due  to  the  change  of  the  communica¬ 
tion  cost  constant. 

Figure  2.6  involves  all  the  vital  information  obtained  in  this  evaluation.  It 
represents  the  2SS’s  minimum  expected  system  cost  vs.  the  optimum  threshold 
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position  and  the  communication  cost  constant. 

More  detailed  observation  of  these  data  is  performed  in  Chapter  5  where  the 
systems  are  compared. 
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CHAPTER  3 


Analysis  of  a  Three-Sensor-System  (3SS) 


3.1.  The  Mode!  and  Configuration 

The  environment  and  other  elements  in  modeling  this  system  are  closely  related 
to  those  in  Chapter  2,  the  two-sensor-system.  The  only  difference  in  this  chapter  is 
that  the  host  sensor  communicates  with  two  slave  sensors,  instead  of  one,  when  the 
host  sensor  makes  an  uncertain  decision.  The  system  cost  is  also  defined  the  same 
way  as  in  (2.2.1)  of  Chapter  2.  Thus,  the  expected  system  cost  expression  is  the 
same  as  (2.3.4).  All  the  expressions  of  terms  in  (2.3.4)  are  directly  applied  for  this 
system. 

3.1.1.  Host  Sensor’s  Decision  Boundaries 

The  design  of  thresholds  in  the  host  sensor  in  the  three-sensor-system  is  very 
similarly  done  as  in  the  two-sensor-system  (refer  to  Figure  3.1).  The  thresholds 
divide  the  observation  space  into  three  decision  regions.  No  Target  (0),  No  Decision 
(?),  and  Target  Detected  (1). 

3.1.2.  The  Overall  Process  of  The  Model 

The  host  sensor’s  confident  local  decision,  either  0  or  1,  will  become  the  final 
decision  of  the  system.  In  case  the  decision  of  the  host  sensor  is  dubious  (yns 
between  TL  and  TU,  or  TL3  and  TU3  in  Figure  3.1),  the  host  sensor  will  request 


binary  informations,  Ussi  and  Uss2.  froin  both  of  the  slave  sensors,  which  are  also 
independently  generated  according  to  their  observation,  yssi  and  yss2>  at  the  slave 
sensor  1  (SSI)  and  the  slave  sensor  2  (SS2).  These  communication  process  also 
involves  a  communication  cost  constant  as  in  2SS.  The  illustration  of  the  process  is 
in  Figure  3.2. 

3.2.  Evaluation  of  Error  Caused  by  Team  Strategies 

^eieam  ^as  the  same  probabilistic  expression  as  that  of  2SS  except  now 

the  expression  is  conditioned  on  two  slave  sensors,  not  one.  The  expression  is  shown 
below. 

Pei^^  =  Pr(crror  resulted  by  communication  using  team  strategy)  =  Pr(E) 


Figure  3.2  Model  Configuration  of  3SS 


=  P,(E  I  yHs  e  [TL.TU])  -P^CyHs  e  [TL.TU]) 

=  P,(E.  y„s  6  [TL.TU]  I  Ho)  ■  P,(Ho)  +  P,(E,  y„s  6  [TL.TU]  I  H, )  ■  Pr(Hi ) 

=  P,(E.  y„s€  [TL.TU]  I  Ho,  Uss,.Uss2)-Pr(Ussi  '  Ho)-P,(Uss2  I  Ho)-P,(Ho) 


29 


+  Pr(E.  yHse  [TL.TU]  1  H,.Ussi.Uss2)-P,(Ussi  I  H,)-P,aJss2  '  Hi)-Pr(H,)  (3.2.1) 


3J.  The  Likelihood  Ratio  Test 

In  evaluating  the  Likelihood  Ratio  Test  (LRT),  it  is  assumed  that  the  individual 
observations  received  at  each  sensor  are  independent  from  each  other.  Thus,  their 
local  decisions  are  also  independent  from  other  sensors.  In  other  words,  the  local 
decision  of  the  slave  sensors  is  statistically  independent  from  (or  not  coupled  with) 
either  the  local  decision  or  the  observation  of  the  host  sensor. 


A(yHS*Ussi*Uss2)  = 


PfCyns.  Ussi.  Uss2  I  Hi) 
Pr(yHS.  Ussi.  Uss2  >  Ho) 


Up=l 

> 

< 

Uf=0 


Using  the  above  assumption,  (3.3.1)  can  be  written  as  following. 


Pr(yHS  I  Hi)-P,(Ussi  I  Hi)-P,(Uss2  ■  Hi)-P,(Hi) 
Pr(yHS  I  Ho)-Pr(Ussi  I  Hq)  ’  Pr(Uss2  i  Ho)-Pr(Ho) 


Uf=1 

> 

< 

Uf=0 


(3.3.1) 


(3.3.2) 


Uf=1 

Pr(yHS  I  Hi)  >  ^  Pr(Ho)  P/UgSl  <  Hp)  P,(Uss2  '  Hp) 
Pr(yHS  I  Hp)  P,(H,)  ■  P,(Ussi  I  H,)  ’  Pr(Uss2  '  H,) 


f(Ussi.Uss2)  =  Tp ■ 


Pr(UsSl  I  Hp) 
Pr(Ussi  I  H,) 


Pr(UsS2  I  Hp) 
Pr(UsS2  I  H,) 


(3.3.3) 


Uf=1 

g(yHs)  <  f(Ussi,Uss2)  (3.3.4) 

Uf=0 


To  is  a  ratio  of  a  priori  probabilities.  The  function  f(Ussi.Uss2)  represents  the  final 
threshold  in  the  host  sensor  after  communication  between  the  host  sensor  and  the 
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slave  sen‘?ors.  As  it  is  seen  in  (3.3.3),  the  function  f(UssitUss2)  have  four  dif¬ 
ferent  values  (thresholds)  depending  on  the  decision  of  the  slave  sensors, 
Ussi  andUss2»  since  the  decisions  of  SS,  Ussi»  Uss2  always  either  0  or  1  . 
This  gives  more  explicit  expression  of  the  function  f(Ussi.Uss2)  which  is  listed 
below. 


.  ,  P,(Ussi=OIHo) 

f(Ussi=0.Uss2-0)-Xt  To  p^(Ussi=OIHi) 

Pr(UsS2=OIHo) 

Pr(Uss2=0IHi) 

(3.3.5) 

,  ^  Pr(Ussi=OIHo) 

f(Ussi-0.Uss2-l)-^t  0  p^(Uss,=OIHi) 

Pr(Uss2=llHo) 
Pr(Uss2=l  IH,) 

(3.3.6) 

. . .  P,(Ussi=llHo) 

f(Ussi-l.Uss2-0)-^t  To  p^(Ussi=llHi) 

Pr(UsS2=4)  1  Ho) 

Pr(UsS2=0IHi) 

(3.3.7) 

. . Pr(UsSl=llHo) 

f(Ussi  =  l.Uss2-l)-^t  To  p^(Ussi=liHi) 

Pr(UsS2=llHo) 

Pr(Uss2=llHi) 

(3.3.8) 

Using  (3.3.4),  we  can  write  the  final  decision.  Up,  rules  of  the  system  as 


Uf  = 


1 

1 

0 

0 


if  yns  ^  TU 

if  g(yHs)  -  f(Ussi>UsS2) 

if  g(yHS)  <  f(UsSl.UsS2) 

if  yHs  <  tl 


Then,  it  is  possible  to  express  Pexjjn,  shown  below. 


P,(E.  y„s  e  ITL.  TU]) 

=  Pr(Ho)-  I  I  Pr(E.  y„se[TL.TU]IUss..Uss2. 

Cs3.=0  l-‘5n=0 
L.'ssi  =  l  Css3=I 


-rP^(Ho)-  2^  X  Pr(E,  Yus  £  (TL, TU)  I  Ussi.  Uss2 


Ho) -P, (Ussi  I  Hq) -PrCUss;  I  Ho) 

,H,)-P,(Uss2  I  H,)-Pr(Uss2  I  H,) 
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Uai»l  Usa“l  Uc-l 

:P.(Ho)-  S  Z  P.lg(yHs)  7 


Un,^  Uatt^O 


f(Ussi»Uss2)  yws  s  [TL,  TU]  I  Ussi. Uss2.l'^)  ‘Pi(Ussi  I  Hq)  •P,(Uss2  I 

Uiji«l  Ua2*l  Uc=0 

■Pr(Hl)-  2  2  Pr(g(yHs)  > 

Ubi-0  Um-O 

f(Ussi.Uss2)andyHse  [TL.TU]  I  Uss,.Uss2.H,) -P^CUssi  I  Hi)-P,(Uss2  1  HO  (3.3.9) 


3.4.  Calculation  of  C  for  Gaussian  Models 


As  in  Chapter  2,  the  Gaussian  distribution  function  is  used  to  give  examples  in 
expressing  C  so  that  it  can  be  evaluated  numerically. 


3.4.1.  Gaussian  Probability  Density  Function 


The  probability  density  functions  are  shown  below.  They  show  PDF  of  "0"  and 
T"  at  HS,SSl.andSS2. 


-(yHS-tJHSo)^ 

fHSo(yHs)=W - e  ^^0 

'^2rtaHSo 

fHS  1  (yus )  =  -1^ - e 

^2kchSi 

~(y«i~Mssio)^ 

fssio(yssi)=  j"  - - e 

^21105510 

“<yssi-Mssii)‘ 

c  ^  s  1  -  2assi, 

fssii(yssi  )  =  -!== - e 

^2rtassii 

-<yss}-^ss:o)^ 

fss2o  (yss2)  =  ■7=—'" - e 

nilOsSTo 


(3.4.1. I) 


(3.4. 1.2) 


(3.4.1.3) 


(3.4. 1.4) 


(3.4. 1.5) 
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fsS2i  (yss2)  = 


>/2jr< 


-<ysst-^ss2i  f 
^*3sszi 


(3.4. 1.6) 


OSS2i 


3.4.2.  Decision  Boundary  of  SSI  &  SS2 


An  equation  of  the  decision  boundary  for  the  SSs  was  derived,  (2.5. 3.5),  in 
Chapter  2.  The  analytical  expression  of  threshold  in  SSI  and  SS2  follows  the  same 
as  in  Chapter  2. 


Hssio  +dssit 
Tssi  = - - + 


o2 


|4SSli  -|4sSlo 


I 

•‘og=|^tss,  • 


Pr(Ho)l 

Pr(H,) 


(3.4.2. 1) 


|4SS2o  +|4SS2i 

Tss2  = - 


(4SS2i  -JiSS2o 


•logej  ^tss 


Pr(Ho) 
^  Pr(Hi) 


(3.4.2.2) 


Then,  the  decisions  at  the  slave  sensors  are  stated  as  below: 


For  SSI, 


and  for  SS2, 


Ussi  = 


0  -  ‘fyssi<Tssi 
1  ,  ifyssi^Tssi 


(3.4.2.3) 


•f  ySS2<TsS2 
ifyss2^Tss2 


(3.4.2.4) 


Since  the  probability  of  error  incurred  by  communication  and  the  final  decision 
boundaries  are  different  from  Chapter  2,  they  are  written  in  this  section.  The  Gaus¬ 
sian  expression  of  the  FTs  are  written: 
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1-0 


f(Ussi=O.Uss2=0)=XfTo-- 


Tssi  -lissio 


''oSlO 


1-Q 


Tssi-lissii 


I-Q 

f(Ussi=0.  UsS2-l)  =  ■  ^0  - - 7 


CfSSh 


Tssi  ~4ssio 


<^SSlo 


Tssi  -MSSli 


®SSli 


1-Q 


Tss2-1^SS2o 


I 


CsS2o 


1-Q 


Tssi-1Ass2i 


I-Q 


OsS2i 


Tss2-liSS2o 


‘jSS2o 


Tssi-^SS2i 


OSS2i 


f(Ussi=I.Uss2=0)  =  A.fTo-' 


Tssi  -4ssio 


<^SSIo 


l-Q 


Tssi  -Mssu 


CfSSli 


’l'ss2  -  M-SS2o 


1 


<JSS2o 


1-Q 


Tssl-^lSS2l 


OSS2i 


f(Ussi=l.  Uss2=l)-^t  'To  ■ 


Tssi-l^SSio 


«^ssio 


TsS2-4SS2o 


^SS2o 


'i - 


Tssi-liSSli 


Ossii 


Tssi  -|J'SS2i 


'^5521 


The  corresponding  probability  of  error  in  the  team  process  is,  then. 


Pey^  =  P, (error  resulted  by  communicating  with  all  SSs  using  team  strategy) 


I  tjii  Tio  „ 

=  F’r(Ho)'j  J  fsSIo(ySSl)dySSI  ■  J  fsS2o(ySS2)dySS2  ■  J  fHSo(yHs)dyns 


Tssi  •• 

/  fssio*^yssi)‘iyssi  ■  /  i^ss2o(yss2)fiyss2  ■  J  OisoCyiisMyiis 

—  Tjij  f(t.‘sjiaO,  Ussa=l) 


(3.4.2.5) 


(3.4.2.6) 


(3.4.2.7) 


(3 .4,2.8) 
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••  Tua  •• 

+  f  fss\o0ss\)<iyss\  •  f  fsS2o(ySS2)dySS2'  {  fHSo(yHs)dyHS 

Ti,,  —  f(U»,»l.Uaa=0) 


+  J  fssio(yssi)dyssi  •  J  fss2o(yss2)dyss2'  J  fHSo(yHs)dyHsr 

Ti,  Ta:  «Ua,.l.Uaj.l) 


Tai  Tin 

I  fssii(yssi)<iyssi  ■  J  fss2i  (yss2)<iyss2  ■ 


J  fHSi(yHs)dyHS 

f(Ua.=0.UiM=0) 


Tai  ••  “ 

+  J  fssii  (yssi)dyssi  ■  J  fss2i(yss2)dyss2'  J  fnsi  (yHs)dyHs 

—  Ta,  f(Uni=0.UB2=l) 


Tm  •• 

I  fssii(yssi)<iyssi  •  /  fss2i(yss2)dyss2’  J  fHSi  (yHs)dyHs 

Tai  Ua3*0) 


f  fssii(yssi)dyssr  J  fss2i  (yss2)dyss2  ■  J  ^hsi  (yHs)<iyHsf 

Tai  Tju  f(t’ai=l.  Uji3=0 


=  Pr(Ho)- 


• 

r 

l-Q 

Tssi  ~USSlo 

•^SSlo 

. 

«■ 

•  1-Q 

[  TsS2-^SS2o 

► .  n 

f(Ussi-0-  Uss2=^*)  “  4hso 

1  <ysS2o 

1 

‘Jhso 

k 

c  J 

»  d 

1-Q 


Tssi  “Ussio 


^SSIo 


Tss2-liSS2o 

•  n 

f(Ussi-0,  Uss2-1)~4hSo 

^^5520 

k  < 

V 

CfUSo 

.  * 

•  1 

r 

' 

- 

Tssi  -l^ssio 

1-Q 

Tssi  -45520 

■Q 

<^ssio 

<^5520 

. 

.  k 

f(Ussi-l.  Uss2-0)~4hSo 
CfHSo 
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•  ^ 

.  o 

*  ' 

TsS2“USS2o 

.A 

f(Ussi  =  l.UsS2-l)“UHSo 

J 

V 

<^SS2o 

V 

<JhSo 

k  J. 

+  Pr(Hi)- 


» 

•  • 

• 

• 

« 

• 

•  • 

1-Q 

Tssi  -  48511 

1-Q 

TsS2-4SS2i 

.  .. 

1-Q 

f(Ussi=^.  Uss2=0)“4hSi 

OSSli 

«ySS2i 

<yHSi 

, 

. 

. 

^  J 

. 

J 

• 

•  • 

* 

* 

•  1 

' 

1-Q 

Tssi  -4ssii 

k  .  A 

TsS2“4SS22 

1-Q 

f(Ussi=0,  Uss2=1)-4hSi 

OSSli 

<JSS2i 

«JHSi 

>  < 

>• 

' 

• 

'  'I 

Tssi-4ssii 

1-Q 

Tssi  -4SS2t 

1-Q 

f(Ussi=l,  Uss2=0)“4hSo 

OSSli 

OSS2i 

“HS,  J 

V 

^  J 

^  J 

, 

c  J 

+  Q 


Tssi  -^SSli 


®SSli 


T'ss2~M^SS2i 


‘JSS2i 


1-Q 


f(Ussi=l,UsS2=l)~|iHSi 


®HSi 


(3.4.2.9) 


Thus,  C  uf  3SS  in  the  Gaussian  case  is  obtained  by  substituting  (3.4.2.9),  (2.5.4. 1), 


and  (2. 5. 4. 2)  into  (2.3.4).  A  program  listing  is  attached  in  Appendix  B. 


3.4.3.  Numerical  Evaluation  of  C 

The  method  of  performing  the  numerical  evaluation  for  C  is  quite  similar  to 
those  done  in  Chapter  2.  The  parameter  values  used  in  this  section  are  as  follows;  All 
the  standard  deviations,  a,  are  set  to  1.0,  the  mean  of  "0"  observation  is  -1.0,  and  the 
mean  of  T"  observation  is  1.0.  The  communication  cost  constant  is  varied  from  0.0 
to  1.0  in  steps  of  0.05.  The  thresholds  in  the  host  sensor  are  moved  away  from  the 
origin.  TL  moves  to  the  negative  direction  and  TU  moves  to  the  positive  direction 
with  the  relationship  of  TU  =  -TL. 


36 


3.4.4.  Comments  on  .Numerical  Evaluation 

The  results  of  the  numerical  evaluation  are  plotted  in  Figure  3.3,  Figure  3.4  and 
Figure  3.5.  Figure  3.3  shows  that  the  curves  of  the  system  expected  cost  over  the 
threshold  positions  with  different  communication  cost  constant  (CCCl).  When 
CCCl  is  greater  or  equal  to  0.55,  the  curves  are  monotonically  increasing,  giving 
minima  at  the  threshold  position  of  0.0.  Figure  3.4  is  an  enlarged  version  of  Figure 
3.3  which  shows  the  minima  of  the  curves  with  CCCl  less  than  0.55  clearly.  Figure 
3.5  can  be  interpreted  that  the  minimum  expected  system  cost  increases  as  CCCl 
increases;  however,  near  CCCl  =  0.5,  the  minimum  expected  system  cost  tends  to  be 
flattening  since  the  thresholds  (TL  and  TU)  are  collapsed  into  the  threshold  position 
of  0.0.  This  is  shown  in  Figure  3.6  where  we  plot  the  minimum  expected  system  cost 
vs.  communication  cost  constant  and  vs.  the  optimal  threshold  position.  More  discus¬ 
sion  of  the  results  are  carried  in  Chapter  5. 
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ccc2;  Two-Sensor  Conununicaiion  Cost  Constant 


HS  Threshold  Position 


Figure  3  J  Expected  System  Cost  vs.  HS  Threshold  Position 


CHAPTER  4 


Analysis  of  a  Two/Three-Sensor-System  (2/3SS) 


4.1.  The  Model  and  Configuration 

The  difference  between  this  chapter  and  the  previous  chapters  is  that  here  the 
host  sensor  chooses  a  communication  scheme,  based  on  quality  of  its  own  decision. 
For  example,  when  the  host  sensor’s  observation,  yns.  falls  in  a  cenain  region  of 
uncertainty,  it  communicates  with  only  one  slave  sensor.  It  communicates  with  two 
slave  sensors  when  the  observation  falls  in  the  other  region  of  uncenainty.  Contrary 
to  the  previous  chapters,  two  different  communication  cost  constants  are  considered; 
one  for  communicating  with  one  slave  sensor,  and  another  for  communicating  with 
two  slave  sensors. 

4.1.1.  The  Host  Sensor’s  Thresholds  and  Decision  Regions 

In  the  host  sensor’s  observation  space,  there  are  four  thresholds  that  divide  the 
space  into  four  decision  regions.  (Actually,  there  are  five  decision  regions  but  two 
out  of  five  regions  yield  the  same  decision.)  When  yn^,  falls  below  TLl  (TL31  in 
Figure  4.1)  or  above  TUI  (TU31  in  Figure  4.1),  the  host  sensor  decides  0  or  1, 
respectively.  When  the  observation  falls  between  TLl  and  TL2  (TL31  and  TL32  in 
Figure  4.1)  or  between  TU2  and  TUI  (TU32  and  TU31  in  Figure  4.1),  the  host 
sensor’s  decision,  Uhs-  becomes  uncenain  (?1).  In  case  of  the  observation  lies 
between  TL2  and  TU2  (TL32  and  TU32  in  Figure  4.1),  the  host  sensor  makes  a  dubi- 
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ous  decision  (?2),  Thus  there  are  level  of  confidence  in  making  uncertain  decision. 
In  ?1  decision  region,  the  probability  of  making  a  correct  decision  is  much  greater 
than  the  probability  of  making  a  false  decision;  then,  a  minimum  help  from  the  slave 
sensors  is  needed.  In  ?2  decision  region,  the  probability  of  making  a  correct  decision 
and  the  probability  of  making  a  false  detection  are  compatible;  thus,  this  situation 
requires  more  information  to  make  a  correct  decision. 

When  the  host  sensor  makes  a  binary  decision  (either  0  or  1),  it  becomes  the 
final  decision  of  the  system.  In  case  the  decision  of  the  host  sensor  is  ?1,  the  host 
sensor  requests  information  only  from  one  of  the  slave  sensors,  say  slave  sensor  1 
(SSI).  On  the  other  hand,  when  the  the  host  sensor  determines  ?2,  it  asks  an  assis¬ 
tance  from  both  of  the  slave  sensors,  slave  sensor  1  (SSI)  and  slave  sensor  2  (SS2). 


Uhs  =  0 

Uhs  =  ?1 

Uhs  = 

j 

Uhs  =  ?1 

1 

TL 

^32 

1 

TL31  TU31 


Figure  4.1  Decision  Boundaries  of  HS  for  2/3SS 
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This  process  is  illustrated  in  Figure  4.2. 


Figure  4.2  Model  Configuration  of  2/3SS 
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4J.  Definition  of  the  System  Cost  Function 

The  following  is  the  cost  function  of  the  system,  C(„„),  which  is  represented  by 
error  probabilities  of  the  individual  sensors.  For  descriptions  of  symbols  used  in  this 
chapter,  please  refer  to  the  beginning  of  this  paper  under  "Symbols  used  in  Chapter 
4". 

C(f)  =  C(zxi  ,Zx2tTLl,TL2,TU2,TUl,Cxi,Cx2) 

=  (1  -  zti  )  •  (1  -  Zt2)  •  Pchs 

+  Zt1 '(1-2x2) '(Pcti  +CTi) 

+  ZX2  •  (1-2x1  )'(PeT2 +CX2)  (4.2.1) 

As  shown  in  the  above  equation,  the  communication  schemes  are  dependent 
upon  the  values  of  zji  and  zj2.  zti  and  z-^  take  binary  numbers  depending  on  the 
type  of  host  sensor’s  uncenain  decision.  When  Uhs  =  ?1.  zji  becomes  1.  When 
Uhs  =  zt2  becomes  1.  This  is  shown  in  the  table  below. 


Communication  Scheme 

Zti 

Zt2 

No  Communication 

0 

0 

Communication  with  SS 1 

1 

0 

Communication  with  SSI  &  SS2 

0 

1 

Table  4.1  Communication  Scheme  of  2/3SS 


4j.  Evaluation  of  an  Expected  System’s  Total  Cost,  C 
Let’s  evaluate  the  expected  cost  of  the  system. 

C  =  E(C(0} 
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—  •Pr(Z7i=0)-Pr(z-r2=0) 

+  C(Olzn»l.  ZT2=0  ■  Pr(ZTl =1)  ■  Pr(ZT2=^) 

+  C(f)lzn-l.*n-I  •Pr(ZTl=l)-Pr(ZT7=l) 

=  PeHs  •P,(ZTl=0)-P,(Z-i7=0) 

+  (Pen  +Cti)  •  Pr(ZTi  =  l)  •  Pr(ZT2^) 

+  (Pen  +Ct2)  •  Pf(ZTl=0)  •  Pr(ZT2=l) 

“Pchs  ■‘■(Pcti  -PeHs)‘Pr(ZTi=l)  +  (Pen'*'‘^T2-PeHs)‘Pr(ZT2=l) 

+  (PeHs  “Pe-n  “Pen  -c-R)'Pr(ZTi=l)-Pr(ZT2=l)  (4.3.1) 

The  terms,  Pens*  Pr(2Ti=l)i  and  Pr(zT2=l).  are  written  in  generalized  probabilistic 
expressions: 

Pbhs  =  Pr (false  local  decision  at  HS) 

=  P, (Decide  H,  I  Ho)  ■P,(Ho)  +  P, (Decide  Hq  I  Hi)-P,(Hi) 

=  P,(y„s  2TU1  I  Ho)-P,(Ho)  +  P,(yHs  ^TLl  I  H,)-P,(H,)  (4.3.2) 

P,(zTi=l)  =  Pr(uncertain  decision  =  ?1;  communication  channel  open  only  with  SSI) 

=  P,(TL1  <  yHs  <  TL2)  +  P,(TL2  <  y„s  TUI) 

=  P,(Ho)  •  {P,(TL1  < y„s  <TL2  I  Ho) •  P,(Ho)  +  P,(TL2  < y„s  <TU1  I  Hq)) 

+  P,(Hi)-(Pr(TLl<y„s<TL2IH,)-P,(Ho)  +  P,(TL2<yHs<TUl  I  H,))  (4.3.3) 

P,(zn=l)  =  Pr(uncertain  decision  =  72;  communication  channel  open  with  SSI  &  SS2) 

=  P,(TL2<yHs<TU2) 

=  P,(TL2  <  yHs  <  TU2  I  Ho)  ■  P,(Ho)  +  P,(TL2  <  y„s  <  TU2  I  H, )  •  P,(H, )  (4.3.4) 
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There  are  two  different  costs  (probability  of  error)  incurred  in  communication, 
Pe-j.j  and  ,  since  the  system  has  two  different  modes  of  communication,  com¬ 
municating  with  one  slave  sensor  (SSI),  or  with  two  slave  sensors  (SSI  and  SS2), 
respectively. 

=P,(error  resulting  after  communication  with  SSI  only)  =  Pr(El) 

=P,(E1  I  VHse  (TLl.TL2]or[TU2.TUl])-P,(yHse  [TLl.TL2]or[TU2.TUl]) 

=  P,(E1.  YHse  [TLl.TL2]orlTU2.TUl]  I  Ho)*P,(Ho) 

+  Pr(El.yHse  [TLl,TL2]or[TU2.TUl]  1  Hi)-P,(H,) 

=P,(E1,  yHs6  n^l.TL2]or(TU2.TUl]  1  Ho,  Ussi)  ■  ?r(Ussi  I  Ho)-Pr(Ho) 

+  P,(El,yHs6[TLl,TL2]orrrU2.TUl]IH,,Ussi)-Pr(Ussi  lH,)-Pr(Hi)  (4.3.5) 

p.^  =  p,(Error  resulting  after  communication  with  SSI  &  SS2)  =  P:(E2) 

=  P,CE2  1  yHs  €  (TL2.'nj2])  ■  P,(yHs  €  [TL2,TU2]) 

=  Pr(E2,  yns  e  [TL2.TU2]  I  Ho)-Pr(Ho)  +  Pr(E2,  yns  e  [TL2.TU2]  I  Hi)-P,(Hi) 

=  P,CE2.  yHs€  [TL2,TU2]  I  Ho.Ussi.UsS2)-Pr(Us3I  '■  Ho)-Pr(Uss2  I  Hor  ?r(.Ho) 

+  P,(E2.  /Hse  [TL2,TLf2]  I  Hi ,  Ujsi .  Uss2)  ’  Pr (Ussi  I  H,)-P,(Uss2  '  Hi)-P,(Hi)  (4.3.6) 


4.4.  The  Likelihood  Ratio  Test 

In  this  chapter,  it  is  necessary  to  evaluate  two  kinds  of  LRT,  since  the  LRT  for 
the  different  communication  schemes  differs.  These  evaluations  closely  follow  those 
derived  in  Chapters  2  and  3. 
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4.4.1.  LRT  for  Communicating  with  SSI  Only 


Let  LR  of  this  case  be 


Ati  (yHS.Ussi)  = 


Pr(yHS.  Ussi  I  Hi) 
Pr(yHS.  Ussi  I  Ho) 


(4.4. 1.1) 


Since  the  observations  received  at  different  sensors  are  mutually  independent, 
(4.4. 1.1)  can  be  written  as 


Ati  (yHs.Ussi) 


Pr(yHS  IHt)-P,(Ussi  IHi)-P,(Hi) 
Pr(yHS  I  Ho)-Pr(Ussi  I  Ho)-P,(Ho) 


Uf=1 

> 

< 

Uf=0 


(4.4. 1.2) 


Thus, 


Pr(yHS  I  Hi)  .  Pr(Ho)  Pr(UsSl  ■  Hp) 

Pr(yHslHo)  ,  ‘  P,(H,)  ^.(Ussi  I  Hi) 


Recalling  the  definitions  made  in  Chapter  2,  (2.4.5),  (2.4.6),  and  (2.4.7),  then 
(4.4. 1.3)  can  be  written  as  below,  provided  that  we  substitute  g(yHs)  =  gTi(yHs)  and 
f(Uss)  =  f(Ussi)- 


Uf=1 

gTi(yHs)  <  f(Ussi)  (4.4. 1.4) 

Uf=0 


The  function  f(Ussi)  represents  the  final  threshold  at  the  host  sensor  after  communi¬ 
cation  with  one  slave  sensor,  SSI.  ffUssi)  can  be  two  different  values  (thresholds) 
depending  on  the  decision  of  the  slave  sensor,  Ussi-  More  explicit  expression  of  the 
function  f(Ussi)  is  listed  below. 
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P,(Ussi=OIHo) 

f(Ussi  -0)  -  Xt  ■  To  •  ^  I  ^ 


(4.4. 1.5) 


Pr(UsSl=llHo) 

f(Ussi  - 1)  =  A-t  ■  To  •  J 


(4.4. 1.6) 


Then,  the  probability  of  error  caused  by  the  team  process  with  SSI  only  can  be 
expressed  in  probabilistic  terms  as  below. . 

Pr(El.  yHse  [TL1.TL2]  or[TU2.TUl]) 

=  Pr(Ho)  •  2  PrOEl.  yHs  6  ITLl.  TL2]  or  [TU2.  TUI]  1  Ussi .  Ho)  •  PrCUssi  I  Ho) 

+  Pr(Ho)-  2  P,(E1,  ynse  [TL1.TL2]  or[TU2.TUl]  I  Ussi.Hi)-Pr(Ussi  I  Hi) 

Ua,-0 

Uai»t 

=  P,(Ho)'  X  Pr{g(yHs)> 

Ua,-0 

f(Ussi)andyHs6  [TL1,TL2]  or  (TU2,TU1]  I  Ussi.Ho) -PtCUssi  I  Ho) 

Cjji=l 

+  Pr(H,)-  X  Pr(g(yus)> 

Ujj,=0 

f(Ussi)andyHsG[TLl.TL2]or[TU2.TUl]  lUssi.Hil-PrlUssi  I  H,)  (4.4.1.7) 


4.4.2.  LRT  for  Communicating  with  SSI  and  SS2 

This  section  is  very  similar  to  the  section  3.3  of  Chapter  3.  The  LRT  when  com¬ 
municating  with  two  slave  sensors  had  been  derived  in  Chapter  3.  Adapting  (3.3.1) 
and  (3.3.2),  we  obtain  AT2(yHS>Ussi»Uss2)  =  •^(yHS«Ussi -Uss:)-  The  (3.3.3)  can  be 
directly  applicable  in  this  section  as  well.  By  replacing  g(yHs)  in  (3.3.4)  with 
2T2(yHs)>  i’nve  the  description  of  the  two-helper  LRT  as 
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Uf=1 

g-r7(yHs)  <  f(Ussi.Uss2) 

Uf=0 

From  (4.4. 1.4)  and  (4.4.2. 1),  the  final  decision.  Up,  rule  of  the  system  can  be  written 
as 


UfH 


1 

1 

1 

0 

0 

0 


.  if  yns  ^  TUI 
.  if  gT2(yHs)  ^  f(UsSl*UsS2) 
.  if  gTl(yHs)  -  fCUsSl) 

.  if  gTl(yHs)  <  f(UsSl) 

•  if  gT2(yHs)  <  f(UsSl.UsS2) 
.  if  yns  <  TLI 


An  explicit  expression  of  the  final  threshold,  f(Ussi»Uss2).  is  dependent  upon  the 
decision  on  the  slave  sensors  as  mentioned  in  Chapter  3.  The  explicit  expressions  are 
given  in  section  3.3,  (3.3.5),  (3.3.6),  (3.3.7),  and  (3.3.S).  Then,  it  is  possible  to 
express  the  probability  of  error  caused  by  using  data  from  two  slave  sensors.  This  is 
shown  in  (4. 4.2.2). 

P,(E2.  yuse  [TL2.TU2]) 

=  Pr(Ho)  •  '  PrCE2.  yHs  G  (TL2,  TU2]  I  L'ssi.  Uss2.  Hq)  •  Pr(Ussi  I  Ho)  •  P,(Uss,  I  Ho) 

Cjsi=l  Cjjj=l 

+  Pr(Ho)-  S  I  P,(E2.yHs€  [TL2.TU2]  IUssi.Uss2.Hi)-Pr(Uss2  I  H,)-P.(Uss2  I  H,) 

U]J23l 

=  Pr(Ho)-  I  I  Pr(g(yHs)> 

i(Ussi  .Uss2)  iind  yiis  g  [TL2,  TU2]  I  Ussi  ■  bss2'Ho) '  PrCUssi  I  Ho)  ■  Pr(iis2  I  Hq) 
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Umi”1 

+  Pr(H,)-  I  I  Pr{g(yHs)> 

Usu~OUni>« 

f(Ussi>Uss2)  and  yHs  €  [TL2,TU2]  I  Ussi.  Ussa.Hi) 'PrCUssi  IHi)-Pr(Uss2  1  Hj)  (4.4.2.2) 

4.5.  Calculation  of  C  under  Gaussian  Models 

We  again  assume  the  probability  density  function  on  the  observation  of  the  host 
sensor  and  the  slave  sensors  given  in  Section  3.4.1  of  Chapter  3.  The  decision  boun¬ 
dary  for  the  local  decision  on  the  slave  sensors  is  also  given  in  Section  3.4.2  of 
Chapter  3,  namely  set  the  binary  decision  threshold  at  0. 


^chs  fHSo(yHs)dyHs +Pr(Hi)-  j  fHSi(yHs)dyHs 


TUl-tiHSo  TLl-pHSi 

=  P,(Ho)-Q  — - ^  +P,(H,)il-Q  — - - i- 

<^HSo  I  <^HSi 

L  J  [  L 


(4.5.1) 


Pr(ZT!=3)  =  Pr(Ho)-j  J  feo (yHs)dyHS  +  J  fHSo(yHs)dyHS 


+  P,(Hi)-|  I  fHSi  (yHS)dyHS  +  J  fHSi  (yHs)dyus 

(tLI  TU2  i 


TLI-uhso'I  TL2-unSol  TU2-UhSo 


=  R(Ho)-  Q  - ^  ~Q|  - ^  +Q  - -Q  - C 

'^HSo  [  tJHSo  Onso  CTiiSo 


"PUl  -uhso'i] 


TL1-uhsi  TL2-UHS,  TU2-UhSi  "n-Il-iiHSi  I 

P,(H,)-  Q  — - ^  ~Q  — - ^  +Q  - - L  -Q  - ^  (4.5.2) 

<JHSi  CTHS,  CThs,  OHSi 
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TU2 


TU2 


Pr(ZT7=l)  =  Pr(Ho) '  J^fHSo(yHs)dyHS  +  Pr(Hi) 


(yHs)dyHs 


=  Pr(Ho)- 


TL2-^hSo 


<JHSo 


-Q 


TU2-|iHSo 


OHSo 


f 

r  y 

-  o 

TL2-|iHSi 

-Q 

TU2-HHS, 

V 

*  y 

OHSi 

ohsi 

r 

« 

L  J 

w  J 

4 

=Pr(Ho)-j 


-  Tss, 

j  fHSo(yHs)dyHs  •  J  fssio(yssi)dyssi 

f(Ussi=0)  ^ 


+  J  fHSo(yHs)dyHs  •  J  fssio(yssi)dyssif 

f(l^'ssi=I)  Tssi 


+  Pr(Hi)- 


f(T-'ssi=>0)  Tssi 

J  fHSi (yHs)dyHs •  J  fssii(yssi)dyssi 


f(Lssi“*)  — 

J  fHSi  (yHs)dyHs  ■  J  fssii(yssi)dyssi’ 

Tssi 


J 


=  Pr(Ho) 


Q 

(• 

f(Ussi=0;-.UHSo 

l-Q 

r  ■> 

Tssi  -^SSlo 

<3HSo 

<^SSlo 

r 

. 

> 

f(UsSl  =  l)-)lHSo 

.  o 

Tssi  “KSSlo 

<^HSo 

.. 

V 

^^ssio 

L  J. 

(4.5.3) 
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• 

r 

f 

+  Pr(Hl)- 

1-Q 

f(Ussi=0)-)iHSi 

»  *4 

1-Q 

* 

OHSi 

. 

Tssi-I^SSli 


Ossii 


1-Q 


f(Ussi=l)-|^HSi 


OHSi 


hQ 


Tssi  "M-SSli 


Ossii 


P.^„=Pr(Ho)- 


(4.5.4) 


r 

» 

• 

» 

1-Q 

Tssi  -l^SSlo 

1-Q 

1. 

TsS2-|^SS2o 

^  *n 

f(Ussi=0.  Uss2=0)  -  I^HSo 

« 

<Jssio 

L.  J 

►  .. 

<^SS2o 

t,  ^ 

V 

CfHSo 

.  ^ 

»  • 

' 

1-Q 

Tssi  -4ssio 

Tss2-11SS2o 

.n 

f(l-lssi=0>  Uss2=l)  — UnSo 

<Jssio 

.  -Q 

Oss2o 

V 

CfHSo 

h 

.  < 

. 

Tssi-4ssio^ 

1-Q 

Tssi-liSS2o 

► .  n 

f(Ussi  =  1 .  Uss2=0)  -  I^HSo 

j 

•< 

c^ssao 

'  V 

<^HSo 

*  > 

+  Q 


Tssi-l^SSio 

r  1 

■rss2  -  liSS2o 

■Q 

^  y 

f(Ussi=l'Uss2=l)“4HSo 

<Jssio 

w  > 

■Q 

<JSS2o 

•  > 

^HSo 

+  Pr(Hi)- 


1-Q 


i-Q 


Tssi  -4ssii 


<Jssii 


• 

* 

1-Q 

’l'sS2-l^SS2i 

>  *4 

1-Q 
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4.5.1.  Numerical  Evaluation  of  C 


>1 


The  same  method  is  used  as  in  the  previous  chapters  in  evaluating  C  numeri¬ 
cally.  The  difference  is  that  the  host  sensor  in  this  system  has  4  thresholds,  unlike 
2SS  and  3SS.  The  thresholds  are  varied  with  a  relationship  of  TU31  =  -TL31,  TU32 

=  -TL32,  aiid  TU32  =  yTU31.  This  threshold  relationship  is  selected  arbitrarily. 

For  the  threshold  configuration,  refer  back  to  Figure  4.1.  The  program  written  for 
this  evaluation  is  attached  under  Appendix  C.  As  shown  in  Figure  4.3,  depending  on 
the  communication  cost  constants,  CCCl  and  CCC2,  individual  curves  are  obtained. 
The  dotted  curve  which  is  generated  using  CCCl  =  0.325  and  CCC2  =  0.65  is  the  last 
curve  with  a  minimum  other  than  at  the  threshold  position  of  0.0.  Figure  4.4  is  an 
enlai'ged  version  of  Figure  4.3.  In  Figure  4.3  and  4.4,  the  curves  have  ripples,  unlike 
the  set  of  curves  shown  in  the  previous  chapters.  This  phenomenon  is  induced  from 
the  arbitrary  choice  of  thresholds,  giving  a  suboptimal  threshold  locations,  and  from 
the  changes  of  the  system’s  communication  scheme  from  one  to  another.  Figure  4.5 
shows  the  minimum  expected  system  cost  holds  at  a  constant  beyond  the  communica¬ 
tion  cost  constant  of  0.65.  This  is  because  that  the  optimum  thresholds  at  the  host 
sensor,  TU31,  TU.32,  TL31,  and  TL32,  eventually  become  zero.  This  is  more  clearly 
represented  in  Figure  4.6. 
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Comparison  of  C  of  2SS,  3SS,  and  2/3SS 


5.1.  Comparison  of  C 

In  this  section  numerically-evaluated  expected  system  costs  in  Chapter  2, 
Chapter  3,  and  Chapter  4  are  compared  against  each  other.  The  comparison  are  made 
based  upon  the  data  obtained  using  Gaussian  models  for  the  different  sensor  systems. 

The  system  expected  costs  are  evaluated  over  the  various  threshold  locations  on 
the  host  sensor’s  observation  space  and  different  communication  cost  constant 
incurred  in  communication  between  the  host  sensor  and  the  slave  sensors.  Data  are 
collected  from  the  results  obtained  through  the  C  expressed  in  terms  of  Q(y)- 
functions.  These  informations  are  plotted  and  attached  at  the  end  of  Chapter  2, 
Chapter  3,  and  Chapter  4.  The  summarized  data  are  tabulated  in  Table  5.1,  Table  5.2, 
and  Table  5.3  in  following  sections. 

5.1.1.  Cof  2SS 

Figure  2.3  shows  that  the  total  expected  system  costs  are  evaluated  as  the  thres¬ 
holds,  TL  and  TU,  are  departing  from  the  origin  with  various  communication  cost 
constant,  CCCl.  As  in  Figure  2.3  or  2.4,  some  of  the  curves  have  minima  other  than 
at  the  threshold  position  of  0.0,  some  don’t.  It  is  roughly  seen  thai  a  minimum  of 
curve  occurs  at  the  threshold  position  of  0.0  when  CCCl  >  0.5.  If  CCCl  >  0.5,  the 
communications  between  sensors  are  prohibited  and  the  final  decision  is  made  by  the 
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host  sensor  alone.  The  exact  value  of  the  communication  cost  constant  that  may  not 
give  minimum  (other  than  the  threshold  position  of  0.0)  is  included  between  0.45  and 
0.50  leaning  more  toward  to  0.45.  The  dotted  curve  indicates  that  the  communication 
cost  constant  is  0.45.  As  the  thresholds  move  away  from  0.0,  the  cost  is  increasing 
beyond  the  optimal  threshold  position.  It  begins  to  stop  increasing  near  the  threshold 
position  of  4.0. 

Figure  2.5  is  a  plot  of  extracted  information  from  Figure  2.3.  It  shows  the 
behavior  of  the  minimum  expected  system  cost  due  to  the  change  of  the  communica¬ 
tion  cost  constant.  As  the  communication  cost  constant  becomes  greater,  the 
minimum  expected  system  cost  increases;  however,  the  cost  starts  saturating  at  CCC 
of  0.45.  The  percentage  change  in  the  expected  system  cost  when  the  communication 
cost  constants  are  varied  from  1.0  to  0.0  is  43.7  %.  This  shows  that  the  communica¬ 
tion  cost  constant  takes  a  very  important  role  in  the  system. 

The  relationship  among  the  minimum  expected  cost,  the  optimal  threshold  posi¬ 
tion,  and  the  communication  cost  constant  is  shown  in  Figure  2.6.  The  numericai 
tabulated  data  are  given  in  Table  5.1.  In  Table  5.1,  when  the  minimum  expected  sys¬ 
tem  cost  is  0.1587,  this  means  there  are  no  communications  between  sensors.  This 
number,  thus,  represents  the  cost  of  making  a  final  decision  by  the  host  sensor  only. 

5.1.2.  Cof3SS 

As  in  the  previous  section.  Figure  3.3  shows  that  the  expected  system  cost  vs. 
the  threshold  position  in  the  host  sensor’s  observation  space.  The  dotted  curve  indi- 
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Communication 

Cost 

Constant 

CCC 

Optimum 

Threshold 

Position 

TU 

Minimum  Expected 
System 

Cost 

0.00 

0.700 

0.0891 

0.05 

0.585 

0.1045 

0.10 

0.490 

0.1175 

0.15 

0.405 

0.1283 

0.20 

0.330 

0.1372 

0.25 

0.265 

0.1444 

0.30 

0.200 

0.1500 

0.35 

0.140 

0.1542 

0.40 

0.085 

0.1570 

0.45 

0.035 

0.1584 

0.50 

0.000 

0.1587 

0.55 

0.000 

0.1587 

0.60 

0.000 

0.1587 

0.65 

0.000 

0.1587 

0.70 

0.000 

0.1587 

0.75 

0.000 

0.1587 

0.80 

0.000 

0.1587 

0.85 

0.000 

0.1587 

0.90 

0.000 

0.1587 

0.95 

0.000 

0.1587 

1.00 

0.000 

0.1587 

Table  5.1  Tabulated  Data  of  2SS 

cates  that  the  curves  with  CCC2  >  0.55  have  minima  at  the  threshold  position  of  0.0. 


The  information  in  Figure  3.3  are  summarized  in  Figure  3.5  and  Figure  3.6.  The 


numerical  tabulated  data  of  these  figures  are  li  fted  in  Table  5.2.  From  the  table  there 


is  51.5  %  difference  in  the  e.xpected  system  cost  when  communication  cost  is  varied 


from  1.0  to  0.0. 
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Communication 

Cost 

Constant 

CCC 

Optimum 

Threshold 

Position 

TU 

Minimum  Expected 
System 

Cost 

0.00 

0.785 

U.0773 

0.05 

0.655 

0.0947 

0.10 

0J5O 

0.1092 

0.15 

0.460 

0.1214 

0.20 

0.380 

0.1315 

0.25 

0.310 

0.1400 

0.30 

0.240 

0.1465 

0.35 

0.180 

0.1516 

0.40 

0.125 

0.1553 

0.45 

0.070 

0.1576 

0.50 

0.015 

0.1586 

0.55 

0.000 

0.1587 

0.60 

0.000 

0.1587 

0.65 

0.000 

0.1587 

0.70 

0.000 

0.1587 

0.75 

0.000 

0.1587 

0.80 

0.000 

0.1587 

0.85 

0.000 

0.1587 

0.90 

0.000 

0.1587 

0.95 

0.000 

0.1587 

1.00 

0.000 

0.1587 

Table  5.2  Tabulated  Data  of  3SS 


5.1.3.  C  of  2/3SS 


The  numerical  tabulated  data  of  Figure  4.3,  Figure  4.4,  and  Figure  4.5  is  in 
Table  5.3.  In  Figure  4.3,  it  is  noted  that  the  dotted  curve  occurs  when  CCCl  =  0.325 
and  CCCl  =  0.65.  The  percentage  change  in  the  expected  system  cost  when  CCCl 
changes  from  0.5  to  0.0,  meaning  CCC2  changes  from  1.0  to  0.0,  is  55.3  %.  It  is 
clearly  shown  in  Figure  4.3  that  the  curves  are  leveling  off  near  the  threshold  position 
of  8.0. 
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Communication 

Cost 

Constant 

with  Two  Sensor 

CCC2 

Communication 

Cost 

Constant 
with  One  Sensors 

CCCl 

Optimum 

Inner  Threshold 

Position 

TU32 

Optimum 

Outer  Threshold 

Position 

TU31 

_ 

Minimum  Expected 
System 

Cost 

0.00 

0.000 

0.5150 

1.030 

0.0712 

0.05 

0.025 

0.4425 

0.885 

0.10 

0.050 

0.3850 

0.770 

0.0964 

0.15 

0.075 

0.3375 

0.675 

0.1072 

0.20 

0.010 

0.2925 

0.585 

0.1169 

0.25 

0.125 

0.2550 

0.510 

0.1255 

0.30 

0.150 

0.2175 

0.435 

0.1330 

0.35 

0.175 

0.1825 

0.365 

0.1396 

0.40 

0.200 

0.1525 

0.305 

0.1452 

0.45 

0.225 

0.1200 

0.1498 

0.50 

0.250 

0.0900 

0.1534 

0.55 

0.275 

0.0625 

0.125 

0.1561 

0.60 

0.300 

0.0350 

0.070 

0.1578 

0.65 

0.325 

0.0075 

0.015 

0.0000 

0.000 

0.1587 

0.75 

0.375 

0.80 

0.400 

0.0000 

0.000 

0.85 

0.425 

0.0000 

0.000 

0.1587  1 

0.90 

0.450 

0.0000 

0.000 

0.1587  i 

0.95 

0.475 

0.0000 

0.000 

1.00 

0.500 

0.000 

Table  5.3  Tabulated  Data  of  2/3SS 


5.2.  Comparison  of  Systems 

Since  each  system’s  numerical  evaluation  results  are  collected,  and  5.1.3,  it  is 


possible  to  carry  out  the  performance  comparison  of  these  systems.  Mainly  th^*  sys¬ 
tems’  expected  cost  and  the  optimal  threshold  position  at  different  communication 
cost  constant  are  considered  for  the  comparison.  The  method  used  to  compare  the 
systems  in  this  section  is  that,  first,  the  2SS  is  compared  with  the  rest  of  systems,  3SS 
and  2/3SS.  Secondly,  the  2SS  is  compared  to  2/3SS.  For  the  convenience,  the  com- 
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municadon  cost  constants  of  0.0  and  0.45  are  chosen  to  be  the  bases  of  comparison. 
The  communication  cost  constant  of  0.0  is  selected  since  it  means  that  there  is  no  risk 
in  communication  between  sensors,  in  other  words,  the  communication  between  the 
host  sensor  and  the  slave  sensors  is  encouraged.  The  communication  cost  constant  of 
0.45  are  chosen  because  it  is  the  largest  communication  cost  constant  of  2SS  which 
gives  an  optimal  threshold  position  other  than  0.0. 

Using  the  table  presented  in  the  previous  sections,  at  the  communication  cost 
constant  of  0.0,  the  expected  system  cost  of  3SS  is  13.24  %  less  than  that  of  2SS. 
Comparing  2SS  to  2/3SS,  2/3SS  outperforms  2SS  by  20.10  %  in  the  expected  system 
cost.  In  comparing  with  2/3SS,  the  outer  threshold  location  is  selected  for  the  com¬ 
parison.  2/3SS  has  47.14  %  larger  width  (or  size)  of  the  dubious  decision  region  in 
systems  observation  space.  This  paragraph  is  summarized  in  Table  5.4. 


- 

Type  of  Sensor  System 

CCCI  =  0.0 

CCC2  =  0.0 

2SS 

3SS 

mss 

Improvement  in 

Expected  System  Cost 

O 

o 

13.24  <=Io 

20.10% 

Improvement  in 

Optimal  Threshold  Position 

p 

47.14  % 

CCCl  =  Communication  Cost  Constant  of  communicating  with  one  sensor 
CCC2  =  Communication  Cost  Constant  of  communicating  with  two  sensors 


Table  5.4  Comparison  of  2SS  to  the  Others  with  CCC=0.0 
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In  aspects  of  the  optimal  threshold  position,  2/3SS  has  a  wider  uncertain  deci¬ 
sion  region  than  3SS  by  31.21  %.  With  CCCl  =  CCC2  =  0.0,  2/3SS  performs  about 
7.90  %  better  than  3SS  in  the  expected  system  cost.  This  is  because  2/3SS  requests 
information  from  the  slave  sensors  more  frequent  than  3SS  since  2/3SS  has  a  wider 
uncertain  decision  region.  This  information  are  contained  in  Table  5.5. 

Now  we  consider  system  improvements  in  the  expected  system  cost  and  in 
optimal  threshold  location  with  the  communication  cost  constant  of  0.45  is  con¬ 
sidered.  In  2/3SS  this  communication  cost  constant  is  used  when  the  host  sensor 
communicates  with  two  slave  sensors;  when  the  host  sensor  communicates  with  onlv 
one  sensor,  the  communication  cost  constant  in  this  case  is  a  half  of  the  prior  case, 
0.225.  In  aspects  of  the  expected  system  cost,  the  difference  of  system  cost  between 
2SS  and  3SS  is  0.5  %  in  favour  of  3SS.  For  the  optimal  threshold  location,  3SS  has  a 
wider  uncenainty  region  by  100  %.  In  comparison  of  the  2SS  to  2/3SS,  2/3SS  per¬ 
forms  better  in  the  expected  system  cost  by  5.43  %.  These  are  listed  in  Table  5.6. 


Type  of  Sensor  System 

CCCl  =  0.0 

CCC2  =  0.0 

3SS 

2/3SS 

Improvement  in 

Expected  System  Cost 

o 

o 

7.90  % 

Improvement  in 

Optimal  Threshold  Position 

0.0% 

31.21  % 

1 

CCCl  =  Communicaiion  Cost  Constant  of  conununicaiir.g  with  one  sensor 
CCC2  =  Communication  Cost  Constant  of  communicating  with  two  sensors 


Table  5.5  Comparison  of  2SS  to  2/3SS  with  CCC=0.0 
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Type  of  Sensor  System 

CCCl  =  125 

CCC2  =  4.5 

2SS 

3SS 

2/3SS 

Improvement  in 
Expected  System  Cost 

0.0% 

0.5% 

5.43  % 

Improvement  in 

Optimal  Threshold  Position 

0.0% 

100% 

584.71  % 

CCCl  =  Communicaiion  Cost  Constant  of  communicating  with  one  sensor 
CCC2  =  Communication  Cost  Constant  of  communicating  with  two  sensors 


Table  5.6  Comparison  of  2SS  to  the  Others  with  CCC  ^  0.0 
2/3SS  performs  about  4.95  %  better  than  3SS  in  the  expected  system  cost.  In 
aspects  of  the  optimal  threshold  position,  2/3SS  has  a  wider  uncertain  decision  region 
than  3SS  by  242.86  %.  This  information  is  contained  in  Table  5.7. 

It  is  noted  that  the  width  of  the  threshold  location  is  shrinking  relatively  faster 
for  2SS  and  3SS  than  2/3SS  as  the  communication  cost  constant  increases.  As  far  as 


Type  of  Sensor  System 

CCCl  =  2.25 

CCC2  =  4.5 

3SS 

2/3SS 

Improvement  in 

Expected  System  Cost 

o 

o 

4.95  % 

Improvement  in 

Optimal  Threshold  Position 

o 

o 

242.86  % 

CCC  I  =  Communicaiion  Cost  Constant  of  communicating  with  one  sensor 
CCC2  =  Communicaiion  Cost  Constant  of  communicating  with  two  sensors 

Table  5.7  Comparison  of  2SS  to  2/3SS  with  CCC  ^  0.0 
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the  expected  system  cost  is  concerned,  there  is  not  a  great  difference  as  in  the  posi¬ 
tion  of  optimal  threshold.  Moreover,  at  a  higher  communication  cost  constant  say 
0.45  (refer  to  Table  5.6),  there  are  insignificant  differences  in  the  expected  system 
cost  among  the  systems. 

It  is  interesting  to  observe  the  relationship  between  optimal  thresholds  and 
Pr(UHS  =  ?)  since  the  probability  of  an  observation  landing  in  the  dubious  region  is 
closely  related  to  the  optimal  thresholds  location  in  the  host  sensor.  PrCUns  = '?) 
represents  that  the  probability  of  the  iioai  sensor's  observation  falls  in  the  uncertainty 
region,  TL  <  yns  ^  TU,  inducing  the  host  sensor’s  local  decision  to  be  This 
relationship  's  tabulated  in  Table  5.8,  Table  5.9,  and  Table  5.10.  It  is  obvious, 
without  looking  at  the  tables,  that  Pr(UHs  =  “?)  decreases  as  the  optimum  threshold 
approaches  to  zero.  When  the  tables  are  plotted  (See  Figure  5.1,  5.2,  and  5.3),  a 
linear  relationship  is  found  between  the  optimal  threshold  positions  and  the  probabil- 


Optimal  Threshold 
Position 

Pr(UHS=?) 

0.700 

0.3375 

0.585 

0.2826 

0.490 

0.2369 

0.405 

0.1959 

0,330 

0.1597 

0.265 

0l282 

0.200 

0.0968 

0.140 

0.0678 

0.085 

0.0411 

0.035 

0.0169 

0.000 

0.0000 

Table  5.8  PrCUns  =  3)  given  the  Optimal  Thresholds  of  2SS 
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ity  of  Uhs  =  ?•  The  program  which  evaluates  the  probability  of  communication  with 
given  the  optimal  thresholds  is  attached  in  Appendix  D. 


Optimal  Threshold 
Position 

P,(Uhs  =  ?) 

0.785 

0.3778 

0.655 

0.3161 

0.550 

0.2658 

0.460 

0.2225 

0.380 

0.1838 

0.310 

0.1500 

0.240 

0.1161 

0.180 

0.0871 

0.125 

0.0605 

0.070 

0.0339 

0.015 

0.0073 

0.000 

0.0000 

Table  5.9  PrCUns  =  V  given  the  Optimal  Thresholds  of  3SS 


Optimal  Threshold 
Position 

Pr(UHS=?) 

1.030 

0.4908 

0.885 

0.4245 

0.770 

0.3707 

0.675 

0.3256 

0.585 

0.2826 

0.510 

0.2465 

0.435 

0.2104 

0.365 

0.1766 

0.305 

0.1476 

0.240 

0.1161 

0.180 

0.0871 

0.125 

0.0605 

0.070 

HUSH 

0.015 

0.000 

0.0000  1 

Table  5.10  PjCUhs  given  the  Optimal  Thresholds  of  2/3SS 
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Figure  5.1  Dubious  Decision  Probability  at  HS  vs.  HS  Optimum  Threshold 


Three-Sensor-System 


Figure  5,2  Dubious  Decision  Probability  at  HS  vs.  HS  Optimum  Threshold 
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Figure  5.3  Dubious  Decision  Probab'.iiy  at  HS  vs.  HS  Optimum  Threshold 


CHAPTER  6 


System  Simulations 


6.1.  Simulation  Method 

Simulation  of  the  systems  evaluated  in  Chapter  2,  3,  and  4  are  performed  to 
understand  how  these  systems  behave  in  a  realistic  environment.  The  same  assump¬ 
tions  as  those  made  in  the  beginning  of  this  work  (Refer  Chapter  1)  are  also  used  in 
the  simulation.  One  additional  system  is  simulated  in  addition  to  three  systems  with 
which  we  have  dealt.  This  system  consists  of  one  sensor  that  has  a  single  threshold 
and  no  slave  sensors. 

To  the  signal,  either  -1  or  1.  Gaussian  noise  is  added  at  the  host  sensor  and  the 
slave  sensors.  In  each  system,  different  slave  sensors  receive  independent  observa¬ 
tion.  However,  in  all  systems,  each  host  sensor  receives  the  same  observation  so  that 
the  performance  of  each  system  can  be  compared  easily.  In  the  simulations,  different 
communication  constants  were  used  and  the  number  of  iterations  performed  was 
10,000.  The  iterations  can  be  interpreted  as  the  number  of  observations  taken  by  the 
host  sensor  and  the  slave  sensors.  In  our  Gaussian  random  number  generation  rou¬ 
tine,  10,000  iterations  provide  with  welt  distributed  Gaussian  random  numbers.  The 
outputs  of  the  different  system  are  compared  in  terms  of  the  percentage  of  correct 
detections  (CD),  false  alarms  (FA),  and  target  misses  (TM).  The  total  detection  error 
is,  then,  FA  -(■  TM.  These  are  listed  in  Table  6.1  and  Table  6.2. 
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Simulation  of  Systems 

Type  of  Sensor  Systems 

ISS 

2SS 

CCC 

CD  (%) 

FA  (%) 

TM(%) 

CD(%) 

FA  (%) 

TM(%) 

0.00 

84.17 

7.64 

8.18 

90.38 

4.60 

5.01 

0.05 

83.71 

8.09 

8.82 

89.02 

5.45 

5.52 

mm 

83.90 

7.96 

8.13 

88.76 

5.51 

0.15 

83.94 

8.00 

8.05 

88.85 

5.60 

li^S 

0.20 

83.89 

8.03 

8.07 

88.26 

5.83 

0.25 

83.32 

8.51 

8.16 

87.24 

6.54 

6.21 

0.30 

83.52 

8.31 

8.16 

86.31 

6.88 

6.80 

0.35 

84.73 

7.66 

7.60 

86.73 

6.72 

6.45 

mm\ 

84.02 

8.11 

7.86 

84.95 

7.69 

7.35 

0.45 

84.03 

8.01 

7.95 

84.70 

7.63 

7.66 

0.50 

84.08 

7.98 

7.93 

84.08 

7.98 

7.93 

mm 

84.45 

7.94 

7.60 

84.45 

7.94 

7.60 

83.88 

8.17 

7.94 

83.88 

8.17 

7.94 

0.65 

84.16 

msm 

7.86 

84.16 

7.97 

7.86 

0.70 

84.81 

7.75  7.43 

84.81 

7.75 

7.43 

CCC  =  communicaiicn  cost  constant 


Table  6.1  Simulation  Results  for  ISS  &  2SS 
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1  Simulation  of  Systems 

]  Type  of  Sensor  Systems 

3SS 

2/3SS 

CCC 

CD  (%) 

FA  (%) 

TM(%) 

CD  (%) 

FA  (%) 

TM(%) 

0.00 

90.36 

2.07 

7.56 

89.93 

3.03 

7.03 

0.05 

89.44 

3.02 

7.53 

90.16 

3.46 

6.37 

0.10 

88.95 

3.34 

7.70 

89.72 

3.53 

6.74 

0.15 

89.38 

3.80 

mmm\ 

90.13 

3.87 

5.99 

0.20 

88.76 

4.28 

KSSI 

89.91 

4.15 

0.25 

87.80 

5.23 

6.96 

89.32 

4.35 

6.32 

0.30 

86.88 

5.73 

7.38 

88.57 

5.12 

6.30 

0.35 

87.19 

5.87 

6.93 

89.16 

4.78 

6.05 

0.40 

85.76 

6.86 

7.37 

87.79 

6.58 

0.45 

84.96 

7.28 

7.75 

87.62 

5.73 

6.64 

0.50 

84.43 

7.74 

7.82 

86.60 

6.50 

6.89 

mm 

84.45 

7.94 

7.60 

86.29 

6.77 

6.93 

0.60 

83.88 

8.16 

7.94 

85.34 

7.32 

mm 

84.16 

7.97 

7.86 

84.48 

7.76 

7.75 

ESQII 

84.81 

7.75 

lAZ 

84.81 

7.75 

7.43 

CCC  3  communication  cost  constant 


Table  6.2  Simulation  Results  for  3SS  &  2/3SS 
6.2.  Simulation  Results  and  Discussion 

The  simulation  results  are  quite  reasonable.  In  general,  the  results  show  that 
2/3SS  performs  the  best  and  followed  by  3SS,  2SS,  and  ISS  in  declining  perfor¬ 
mance.  As  shown  in  the  previous  chapter  the  optimal  threshold  of  2SS  collapses  to 
0.0  when  the  communication  cost  constant  is  0.5.  This  is  also  shown  in  Table  6.1  as 
the  CD  of  2SS  equals  to  that  of  ISS  when  CCC  becomes  0.50.  Also,  CD  of  3SS  and 
2/3SS  become  that  of  ISS  when  CCC  is  equal  to  0.55  and  0.70,  respectively.  This  is 
because  the  host  sensors  in  the  different  systems  receive  the  same  observations. 

There  are  about  6  %  difference  in  CD  between  ISS  and  other  systems;  however, 
the  difference  among  the  2SS,  3SS,  and  2/3SS  is  rather  insignificant  when  CCC  is 
0.0.  This  is  because  there  are  no  influence  of  communication  cost  constant  to  each 


system.  As  the  CCC  increases,  the  differences  in  CD  among  systems  become 
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noticeable  even  though  the  largest  difference  are  about  3.5  %.  For  ISS  CD,  FA,  and 
TM  are  virtually  remain  constant  over  all  the  CCCs  were  used  since  the  performance 
of  ISS  is  independent  from  CCC  CD  in  2SS  decreases  as  CCC  increases.  FA  and 
MT  are  increasing  as  CCC  increases.  In  3SS  and  2/3SS  FA  is  about  3.5  and  2.3  times 
less  than  TM,  respectively,  when  CCC  =  0.0.  As  CCC  increases  the  ratio  of  FA  and 
MT  approaches  to  1.0.  These  informations  are  tabulated  in  Table  6.1  and  Table  6.2. 
The  plotted  version  of  these  data  are  in  Figure  6.1,  Figure  6.2,  and  Figure  6.3.  The 
program  for  this  simulation  is  attached  in  Appendix  E. 
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>  ISS;  X— ->2SS;  A — >  3SS;  □ 


CHAPTER  7 


Conclusion 

One  objective  of  this  study  was  to  characterize  team  strategy  decision  methods 
in  terms  of  analytical  derivation,  numerical  evaluations,  and  system  simulations.  The 
optimization  of  the  system  in  terms  of  minimization  of  either  the  expected  system 
cost  or  the  probability  of  error  in  decision  is  another  objective. 

The  team  strategy  is  applied  to  three  different  systems,  and  the  performance  of 
each  system  is  characterized.  Cost  functions  for  each  system  are  defined.  From  the 
cost  function,  the  expected  system  cost,  C,  is  derived.  The  C  is  represented  in  gen¬ 
eral  probabilistic  terms  as  well  as  for  Gaussian  statistics  using  Q(y)  functions.  The 
numerical  evaluations  are  performed  for  Gaussian  models.  The  numerical  evaluation 
shows,  subject  to  communication  cost,  that  2/3SS  is  the  most  efficient.  The  next  most 
desirable  system  is  3SS  and  the  least  is  2SS.  Simulation  of  the  three  systems  was 
also  carried  out.  The  simulation  results  confirm  the  above  order  of  desirability. 

The  communication  cost  constant  plays  an  important  role  in  the  global  decisions 
of  team  strategies.  Clianges  in  ilie  coramunication  cost  constant  influence  the  fre¬ 
quency  of  communication  between  sensors.  Selection  of  optimum  thresholds  in  the 
host  sensor  is  ueavily  dependent  upon  the  communication  cost  constant  since  the  fre¬ 
quency  of  the  communication  allowance  determines  the  optimal  threshold  positions. 

The  simulation  results  show  that  there  is  a  communication  cost  constant  which 
makes  all  the  systems  perform  the  same.  The  communication  cost  constant  in  this 
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situation  is  0.7.  Thus  wnen  the  communication  oetween  sensors  becomes  very  risky 
or  expensive,  meaning  a  high  probability  of  interception,  the  sensors  avoid  communi¬ 
cation.  This  fact  is  shown  by  comparing  the  system  to  ISS  because  ISS  has  no  com¬ 
munication  capability,  i.e,  there  are  no  slave  sensors  involved.  Refer  to  Figure  6.1, 
6.2,  and  6.3. 

As  the  communication  cost  increases,  the  system  with  three  sensors  apparently 
make  better  chances  of  detection  than  the  system  with  two  sensors.  When  the  com¬ 
munication  cost  is  large  so  that  communications  between  sensors  are  prohibited,  then 
the  performance  of  2SS,  3SS,  and  2/3SS  is  compatible  to  that  of  ISS. 
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APPENDIX  A 


Program  for  Cost  Evaluation  of  2-Sensor-System 


This  appendix  contains  a  FORTRAN  program  listing  which  evaluates  (2.3.4) 
substituted  with  (2.5.4. 1),  (2.5.4.2),  and  (2.5.4.3)  numerically.  It  consists  of  a  main 
routine  called  "TWOSENSYS"  and  two  subroutines,  "QfunE"  and  "FindMin".  Pro¬ 
gram  TWOSENSYS  (TWO  SENsor  SYStem)  is  responsible  for  iterations  (generation 
of  Gaussian  observations),  main  calculations,  calling  subroutines,  writing  outputs  to 
files,  etc.  Subroutine  QfunE  evaluates  Q(y)-function  (refer  to  (2.5. 1.1))  when  limits 
of  integration  are  provided.  Subroutine  FindMin  sorts  a  minimum  in  output  data. 
This  routine  is  used  to  find  an  optimal  threshold  in  HS  where  the  expected  system 
cost  is  minimum. 

The  program  computes  expected  system  costs  over  threshold  positions  in  HS  for 
a  given  communication  cost  constant  Descriptions  of  variables  used  and  comments 
are  embedded  in  the  program. 

Figure  2.3,  Figure  2.4,  Figure  2.5,  and  Figure  2.6  are  plotted  version  of  outputs 
from  this  program.  The  taoulated  data  are  contained  in  l  able  5.1. 
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PROGRAM  TWOSENSYS 


^****a»*»»**»*»****»«»««*****ai«*»*»«*«*«****»«*»*»«**«**«*4i*»***«»»*»** 

c  Author:  Howard  C.  Choe 
c  Organization:  Department  of  Electrical  Engineering 
c  The  University  of  Virginia,  Charlottesville 

c  This  program  numerically  evaluates  expected  costs  of  a 
c  two-sensor-system  which  uses  team  strategies  for  Gaussian  statistics. 

c  Variable  Descripuon 

c  pO  :  a  priori  probability  of  "No  target  exists",  HO 
c  pi  :  a  priori  probability  of  "Target  exists”,  HI 

c  For  Host  Sensor 

c  mhO  :  mean  value  of  HO  received  by  HS 
c  mhl  :  mean  value  of  HI  received  by  HS 
c  shO  :  standard  deviation  of  HO  received  by  HS 
c  shl  :  standard  deviation  of  HI  received  by  HS 
c  For  Slave  Sensor 

c  msO  :  mean  value  of  HO  received  by  SS 

c  msl  :  mean  value  of  HI  received  by  SS 

c  ssO  :  standard  deviation  of  HO  received  by  SS 
c  ssl  :  standard  deviation  of  HI  received  by  SS 
c  Thresholds 
c  For  Host  Sensor 
c  TL  :  lower  Jtreshold 

c  TU  :  upper  threshold 

c  For  Slave  Sensor 
c  Tss  :  LRT  optimal  threshold 
c  Pre-cost  constant 
c  cOO  :  deciding  HO  given  HO 
c  clO  :  deciding  HI  given  HO 
c  cll  :  deciding  HI  given  HI 
c  cOl  :  deciding  HO  given  HI 
c  Team  effort  cost 

c  ccc  :  communication  cost  constant 

c  Declaration 


REAL  pO,  pi 

REAL  mhO,  mhl,  msO,  msl 
REAL  shO,  shl,  ssO,  ssl 


REALcOO.clO.cll.cOl 
REAL  Tss 

REAL  ccc(21),  cfmin(21),  opthr(21) 

REAL  TU(lOOl).  TL(lOOl) 

REAL  pehs(lOOl).  pzcom(lOOl).  cf(1001) 

REAL  cbar(21,1001) 

DATA  ccc/O.OO,  0.05,0.10,  0.15,  0.20, 

•  0.25,0.30,0.35,0.40,0.45, 

•  0.50, 0.55, 0.60, 0.65, 0.70, 

•  0.75, 0.80, 0.85, 0.90, 0.95, 1.00/ 

c  The  input  data 

c  a  priori  probability  of  the  binary  hypothesis  environment 
p0  =  0.5 
pi  =0.5 

c  The  statistics  of  the  received  information 
mhO  =  -i;0 
mhl=  1.0 
msO  =  -1.0 
msl  =  1.0 

sigma  =  1.0 

c  Assign  all  standard  deviation  to  the  same  value 
shO  =  sigma 
shl  =  sigma 
ssO  =  sigma 
ssl  =  sigma 

c  pre<ost  values 
cOO  =  0.0 
cl0=  1.0 
cll  =0.0 
cOl  =  1.0 

c  Evaluation  of  the  prc-calculated  threshold  for  the  team  strategy  and 
c  for  the  slave  sensor 

plamss  =  (clO  -  c00)/(c01  -  cll) 
plamt  =  plamss 

c  Evaluation  of  the  ratio  between  a  priori  probabilities  and  LRT  threshold 
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c  for  the  host  sensor  and  the  slave  sensor,  considering  each  sensor  is 
c  centralized  individually. 

tO  =pO/pl 

Ths  =  (mhO  +  mhl)/2.0  +  (sigma**2/(mhl  -  mhO))*LOG(plamt*tO) 
Tss  =  (msO  +  msl)/2.0  +  (sigma**2/(msl  -  msO))*LOG(plamss*lO) 

c  Evaluation  of  the  final  threshold  after  communication 

aOs  =  (Tss  -  msO)/ssO 
als  =  (Tss  -  msl)/ssl 

CALL  QfunE(aOs,  QaOs) 

CALL  QfunE(als,  Qals) 

fusO  =  plamss  *  tO  *  (1.0  -  Qa0s)/(1.0  -  Qals) 
fusl  =  plamss  *  tO  *  QaOs/Qals 

c  Evaluation  of  Q{y)-function  values  with  fusO  and  fusl 

fusOhO  =  (fusO  -  mh0)/sh0 
fuslhO  -  (fusl  -  mh0)/sh0 
fusOhl  =  (fusO  •  mhl)/shl 
fuslhl  =  (fusl  ■  mhl)/shl 

C.ALL  QfunE(fus0h0,  QfusOhO) 

CALL  QfunE(fuslhO.  QfuslhO) 

CALL  QfunE(fusOhl,  QfusOhl) 

CALL  QfunE(fus  1  h  1 ,  Qfus  1  h  I ) 

c  Calculation  of  an  error  probability  by  the  team  effort 

pcteam  =  (Qfus0h0*(1.0-Qa0s)  +  QfuslhO*QaOs)  *  pO 
•  +  ((1.0-Qfus0hl)»(1.0-Qals)  +  (1.0-Qfuslhl)‘Qals)  *  pi 

c  Calculauon  of  ihc  host  sensor  error  probability  and  that  of 
c  communication  probability,  pchs  and  pzcom,  respectively, 
c  These  probabilities  are  TL  and  TU  dependent  which  means  that 
c  whenever  TL  and  TU  changes,  values  of  pchs  and  pzcom  also  change. 


tine  =  0,02 
DO  10ia=  1.21 
ib  =  0 

DO  30  thr  -  Ths,  Ths  +  4.0.  tine 
ib  =  ib  +  1 
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TU(ib)  =  thr 

TL(ib)  =  Ths  -  (FLOAT(ib)  -  1.0)*tinc 

c  Evaluation  for  pehs  for  various  TL  and  TU 

tuhO  =  (TU(ib)  -  mhOVshO 
tlhl  =  (TL(ib)  -  mhiyshl 

CALL  QfunEftuhO,  QtuhO) 

CALL  QfunEftlhl.Qdhl) 

pehs(ib)  =  QtuhO’pO  +  (1.0  -  Qtlhl)*pl 

c  Evaluation  for  pzcom  for  various  TL  and  TU 

tlhO  =  (TL(ib)  -  mhOyshO 
tuhl  =(TU(ib)- mhiyshl 

CALL  QfunE(UhO.  QtlhO) 

CAIl,  QfunE(tuhl,  Qluhl) 

pzcom(ib)  *  (QtlhO  -  Qtuh0)*p0  +  (Qtlhl  -  Qtuhl)*pl 

c  Evaluation  for  the  expected  probability  of  error  of  the  system,  COST 

cbar(ia,ib)  =  pchstib)  +  (peteam-pehs(ib)+ccc(ia))*pzcom(ib) 
cf(ib)  =  cbar(ia,ib) 

30  CONTINUE 

c  Extract  the  minimum  system  cost  for  each  case  of  ccc(ia) 

CALL  nNDMIN(cf.ib.mini) 
cfmin(ia)  =  cf(mini) 
opthrtia)  =  TL'(mini) 

10  CONTINUE 

c  Write  OUTPUT  DATA . 

WRITE  (10,*)  'Rauo  of  a  prion  probabilities:  tO 
WRITE  (10,*)  "The  LRT  threshold  of  Slave  ‘^ensor:’,  Tss 
WRITE  (10,*)  ’Final  Threshold  when  Us  =  0:  ’,  fusO 

WRITE  (10,*)  ’Final  Threshold  when  Us  =  1;  ’,  fusl 

WRITE  (10,*)  ’Error  caused  by  team  strategies:  peicam 
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DO  50  ic  ==  1.  ib 

WRITE  (1 1,1000)  TL(ic),  TU(ic).  (cbar(id.ic).id=l.lO) 
50  CONTINUE 

DO  70  ie=  l.ib 

WRITE  (12,1000)  TL(ie),  TU(ie),  (cbar(ig,ie),ig=l  1,20) 
70  CONTINUE 

DO  90  ih=  1,21 
WRITE  (13,*)  ccc(ih),  cfmin(ih) 

90  CONTINUE 


DO  linii=  1,21 
WRITE  (14,*)  opthr(ii),  cfmin(ii) 
110  CONTINUE 


DO  130  ij  =  1,21 
WRITE  (15,*)  ccc(ij),  opthr(ij) 

130  CONTINUE 

c  Format  staiemenis 

1000  FORMAT  (”,F6.3,1X,11(F6.4,1X)) 

STOP 

ENT) 


SUBROUTINE  QfunE(x.x.  erfcx) 

c  This  funcuon  calculates  the  error  function  and  the  cor.iplimeniary 
c  error  funcuon  for  the  value  "xx" 

c  Accuracy  is  to  within  1.5E-07. 

REAL  X,  XX,  erfcx 

REAL  al,  a2,  a3,  a4,  a5,  p,  pi,  t 

pi  =  3.141592654 
X  =  ABS(xx) 

al  =  0.319381530 
a2  =  -0.356563782 
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a3  =  1.781477937 
a4  = -1.821255978 
a5=  1.330274429 

p  =  0.2316419 
t=  1.07(1.0  + p*x) 

si  =  al*t  +  a2*i**2  +  a3*i**3  +  a4*t**4  +  a5*l**5 
s2  =  sl''EXP(-(x**2)/2.0) 

IF  (XX  .GE.  0.0)  THEN 
erfcx  =  s2/SQRT(2.0*pi) 

ELSE  IF  (XX  .LT.  0.0)  THEN 
erfcx  =  1.0  -  s2/SQRT(2.0*pi) 

END  IF 


RETURN 

END 


SUBROLrrLNEFLNDMIN(array.isize,minindex) 

REAL  aiTay(isize) 

INTEGER  minindex 


int  =  1 

1 1  CONTIMJE 

DO  10  i  =  ini+l,  isize 
IF  (aiTay(int)  .GT.  airayli))  THEN 
minindex  =  i 
int  =  minindex 
GOTO  1 1 
END  IF 

10  CONTINUE 


RETLEN 

ENT) 


APPENDIX  B 


Program  for  Cost  Evaluation  of  3-Sensor-System 


In  Appendix  B,  a  FORTRAN  program  used  for  numerical  evaluations  of  (2.3.4) 
substituted  with  (2.5.4. 1),  (2.5. 4.2),  and  (3.4.2.9)  is  contained.  The  main  program  is 
called  "THREESENSYS"  (THREE  SENsor  SYStem).  As  is  in  Chapter  2,  the  same 
subroutines,  QfunE  and  FindMin,  are  also  used. 

Outputs  of  this  program  are  represented  by  Figure  3.3,  Figure  3.4,  Figure  3.5, 
and  Figure  3.6.  Some  of  the  data  are  also  available  from  Table  5.2. 


PROGRAM  THREESENSYS 


c  Author  Howard  C.  Choe 
c  Organization:  Department  of  Electrical  Engineering 
c  The  University  of  Virginia.  Charlottesville 

c  Purpose:  Master  of  Science  Research 

c  This  program  numerically  evaluates  expected  system  costs  of  a 
c  three-sensor-system  which  uses  team  strategies  for  Gaussian  statistics. 

c  Variable  Description 

c  pO  :  a  priori  probability  of  "No  target  exists",  HO 
c  pi  :  a  priori  probability  of  "Target  exists",  HI 
c  For  Host  Sensor 

c  mhO  :  mean  value  of  HO  received  by  HS 
c  mhl  :  mean  value  of  HI  received  by  HS 
c  shO  :  standard  deviation  of  HO  received  by  HS 
c  shl  :  standard  deviation  of  HI  received  by  HS 
c  For  Slave  Sensor  1 

c  mslO  :  mean  value  of  HO  received  by  SSI 
c  msl  1  :  mean  value  of  HI  received  by  SS 1 
c  sslO  :  standard  deviation  of  HO  received  by  SSI 
c  ssl  1  :  standard  deviation  of  HI  received  by  SSI 

c  For  Slave  Sensor  2 

c  ms20  :  mean  value  of  HO  received  by  SS2 
c  ms2l  :  mean  value  of  HI  received  by  SS2 
c  ss20  :  standard  deviation  of  HO  received  by  SS2 
c  ss21  ;  standard  deviation  of  HI  received  by  SS2 
c  Thresholds 
c  For  Host  Sensor 
c  TLO  :  lower  threshold 

c  TUQ  :  upper  threshold 

c  For  Slave  Sensor  1 
c  Tssl  :  LRT  optimal  threshold 
c  For  Slave  Sensor  2 
c  Tss2  :  LRT  optimal  threshold 
c  Pre-cost  constant 


c 

cOO 

:  deciding  HO  given  HO 

c 

clO 

:  deciding  HI  given  HO 

c 

cll 

;  deciding  HI  given  HI 

c 

cOl 

:  deciding  HO  given  HI 

c  Team  effort  cost 

c  rrrj)  •  ccrmvr:;  anon  cost  consuni  tor  communicating  with  two  sensor 
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c  Other  variables 

c  cfminO  :  minimum  cost  in  a  single  case  of  run,  i.e.,  for  a  cccQ 
c  opihrO  :  threshold  where  cfminQ  is  occured 

c  pehsO  :  expected  error  of  the  system  when  there  is  no  communication 
c  pzcomO  ■-  probability  of  communication  would  occur 
c  cfO  :  expected  cost  of  the  system  at  various  of  thresholds 
c  cbar(,) :  same  as  cfO  but  saved  in  2-D  array 

c  Other  variables  are  commented  as  program  is  progressed. 

g**«»***»»**** «**«******•»«»****«*•»»*»*»•»«•**:*«»*««*»•**********«« ******* 

c  Declaration 
REAL  pO,  pi 

REAL  mhO,  mhl,  mslO,  msll,  ms20,  ms21 
REAL  shO,  shl,  sslO,  ssl  1,  ss20,  ss21 
REAL  cOO.clO,  ell,  cOl 
REAL  Tssl,  Tss2 

REAL  ccc(2l) 

REAL  cfmin(21),  opihr(2I) 

REALTU(lOOl),  TL(lOOl) 

REAL  pehs(lOOl),  pzcom(lOOl).  cf(1001) 

REAL  cbar(21,1001) 

DATA  ccc/O.OO,  0.05, 0.IO,  0.15, 0.20, 

•  0.25,0.30,0.35,0.40,0.45. 

•  0.50,0.55,0.60,0.65,0.70, 

•  0.75, 0.80, 0.85, 0.90, 0.95, 1.00/ 

c  Determination  of  communication  cost  constant  for  communicating 
c  with  2  slave  sensors 

c  The  input  data 

pO  =  O.j 
pi  =0.5 

mhO  =  -1.0 
mhl  =  1.0 
mslO  =  -1.0 
msl  1  =  1.0 
ms20  =  -1.0 
ms21  =  1.0 
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sigma  =  l.C 


c00  =  0.0 
clO=  1.0 
cl  1  =  0.0 
cOl  =  1.0 

c  Assign  all  standard  deviation  to  the  same  value 

shO  =  sigma 
shl  =  sigma 
ssl0  =  sigma 
ssl  1  =  sigma 
ss20  =  sigma 
ss21  =  sigma 

c  Hva!  nation  of  the  pre -calculated  threshold  for  the  team  strategy  and 
c  for  the  slave  sensor 

plamssl  =  (clO  -  c00)/(c0l  -  cll) 
plamss2  =  (c  10  -  c00)/(c0l  -  cl  1) 
plamt  =  plamssl 

c  Evaluation  of  the  ratio  between  a  priori  probabilities  and  LRT  threshold 
c  for  the  host  sensor  and  the  slave  sensor,  considering  each  sensor  is 
c  centralized  individually. 

tO  =  pO/pl 

Ths  =  (mhO  +  mhl)/2.0  +  (sigma**2/(mhl  -  mhO))*LOG(plamt*tO) 
Tssl  =  (mslO -t- msll)/2.0 

•  +  (sigma**2/(msl  1  -  mslO))*LOG(pIamssl*tO) 

Tss2  =  (ms20  +  ms21)/2.0 

•  (sigma**2/(ms21  ms20))*LOG(plamss2*tO) 

c  Slave  sensor  1 

c  bslO  :  integration  limit  for  Q-function  under  HO 
c  bsl  1  :  integration  limit  for  Q-function  under  HI 
c  QbslO  :  probability  of  making  ussl=l  under  HO 
c  Qbsll  :  probability  of  making  ussl=I  under  HI 

bslO  =  (Tssl  -  msl0)/ssl0 
h<:n  =(Tssl  -msin/scll 


CALL  QfunEfbslO,  QbslO) 
CALL  QfunE(bsll,  Qbsll) 
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c  Slave  sensor  2 

c  bs20  :  integration  limit  for  Q-function  under  HO 
c  bs21  ;  integration  limit  for  Q-function  under  HI 
c  Qbs20 :  probability  of  making  uss2=l  under  HO 
c  Qbs21 :  probability  of  making  uss2=I  under  HI 

bs20  =  (Tss2  -  ms20)/ss20 
bs21  =  (Tss2  -  ms2l)/ss21 

CALL  QfunE(bs20,  Qbs20) 

CALL  QfunE(bs2l.Qbs21) 

c  FOR  COMMUNICATING  WITH  2  SLAVE  SENSORS 
c  Evaluation  of  the  final  threshold  after  communication 
c  fOO  ;  final  threshold  when  ussl=0  and  uss2=0 
c  fOl  :  final  threshold  when  ussl=0  and  uss2=l 
c  flO ;  final  threshold  when  ussl=l  and  uss2=0 
c  fl  1 :  final  threshold  when  ussl=l  and  uss2=l 

roO  =  plamssI*t0*(1.0-Qbsl0)*(1.0-Qbs20)/((1.0-Qbsl  1)*(L0-Qbs2l)) 
fOl  =  plamssl*i0*(1.0-Qbsl0)*Qbs20/((1.0-Qbsll)*Qbs21) 
no  =  plamss2*t0*Qbsl0*(1.0-Qbs20)/(Qbsll*(L0-Qbs21)) 
ni  =  plamss2*t0*Qbsl0*Qbs20/(Qbsll*Qbs21) 

c  Evaluation  of  Q(y)-funcuon  values  using  the  above  values 
c  roOhO  :  integration  limit  with  fOO  under  HO 
c  fOlhO  :  integration  limit  with  fOl  under  HO 
c  f  lOhO  ;  integration  limit  with  flO  under  HO 
c  fl  IhO  ;  integration  limit  with  fll  under  HO 

IDOhO  =  (100  -  mhOyshO 
10Ih0  =  (rai  -mhOyshO 
flOhO  =  (flO  -  mhOyshO 
fllhO  =  (fll-mhOyshO 

c  QfOOhO  :  probability  of  making  uf=l  with  fOO  under  HO 
c  QfOlhO  :  probability  of  making  uf=l  with  fOl  under  HO 
c  QflOhO  :  probability  of  making  uf=l  with  flO  under  HO 
c  Qfl  IhO  :  probability  of  making  uf=l  with  fl  1  under  HO 

CALL  QfunE(f00h0,  QfOOhO) 

CALL  QtunE(fOIfiO,  QlOlhO) 

CALL  QtunE(flOhO,  QflOhO) 

CALL  QfunE(fl  IhO,  Qfl  IhO) 
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c  fDOhO  ;  imegraiion  limit  with  fOO  under  HI 
c  fOlhO  :  integration  limit  with  fOl  under  HI 
c  flOhO  ;  integration  limit  with  flO  under  HI 
c  fl  IhO  :  integration  limit  with  fl  1  under  HI 

fDOhl  =  (nX)-mhl)/shl 
fDihl  =  (roi  -mhiyshl 
flOhl  =  (flO-mhl)/shl 
fllhl  =  (fll  -mhlVshl 

c  QfDOhl  :  probability  of  making  uf=l  with  fOO  under  HI 
c  QfOlhl  ;  probability  of  making  uf=l  with  fDI  under  HI 
c  QflOhl  :  probability  of  making  uf=l  with  flO  under  HI 
c  Qfl  Ihl  :  probability  of  making  uf=l  with  fl  1  under  HI 

CALL  QfunE(fDOhl.  QiDOhl) 

CALL  QfunE(lDlhl,  QlDlhl) 

CALL  QfunE(flOhl,  QflOhl) 

CALL  QfunE(fl  Ihl.  Qfl  Ihl) 

c  Calculation  of  an  error  probability  by  the  team  strategy 
c  with  communicating  with  2  slave  sensors 

pctcam  =  (  Qro0h0*(1.0.Qbsl0)*(1.0-Qbs20) 

•  +Qroih0*(l.O-Qbsl0)*Qbs20 

•  +Qfl0h0*Qbsl0*(1.0-Qbs20) 

•  +  Qf  1 1  hO*Qbs  10*Qbs20  )  •  pO 

•  +(  (I.0-Qf00hl)*(1.0-Qbsll)*(1.0-Qbs21) 

•  +(1.0-Qroihl)*(1.0-Qbsll)*Qbs21 

•  +(1.0-Qfl0hl)*Qbsll*(1.0-Qbs21) 

•  +(1.0-Qfllhl)*Qbsll*Qbs21  )*pl 

c  Calculauon  of  the  host  sensor  error  probability,  and 
c  that  of  communication  frequency  probability  with  2  sensors-pehs, 
c  pzcom-respectively. 

tine  =  0.02 
iX)  10ia=  1,21 
ib  =  0 

DO  30  thr  =  Ths,  Ths  +  4.0,  tine 
ib  =  ib  +  I 
TU(ib)  =  thr 

TL(ib)  =  Ths  -  (FLOAT(ib)  -  1.0)*tinc 
c  Evaluation  for  pchs  for  various  TL  and  TU 


87 


luhO  =  (TU{ib)  -  mhOVshO 
tihl  =  (TL(ib)  -  mhiyshl 

CALL  QfunE(tuhO,  QtuhO) 

CALL  QfunE(tlhl.  Qilhl) 

pehs(ib)  =  QtuhO’pO  +  (1.0  -  Qahl)*pl 

c  Evaluation  for  pzcom  for  various  TL  and  TU 

UhO  =  CTL(ib)  -  mhOyshO 
tub  I  =(TU(ib)  -  mhiyshl 

CALL  QfunE(tlhO,  QtlhO) 

CALL  QfunE(tuhl,  Qtuhl) 

pzcom(ib)  =  (QtlhO  -  QtuhO)*pO  +  (Qtlhl  -  Qtuhl)*pl 

c  Evaluation  for  the  expected  probability  of  error  of  the  system,  COST 

cbar(ia.ib)  =  pehs(ib)  +  (petcam-pehs(ib)+ccc(ia))*pzcom(ib) 
cf(ib)  =  cbar(ia.ib) 

30  CONTINUE 

c  Extract  the  minimum  system  cost  for  each  case  of  ccc(ia) 

CALL  FINDMIN(cf.ibjiiini) 
cfmin(ia)  =  cf(mini) 
opthr(ia)  =  TU(mini) 

10  CONTINUE 

c  Write  OUTPUT  DATA . 

WRITE  (10,*)  'Ratio  of  a  priori  probabilities:  ’,  lO 

WRITE  (10,*)  'LRT  Threshold  for  SSI - :  ’,  tssi 

WRITE  (10.*)  'LRT  Threshold  for  SS2 - :  tss2 

WRITE  (10,*)  'FT  at  HS  when  Usl=0  and  Us2=0  fOO 
WRITE  (10,*)  'FT  at  RS  when  Us  1=0  and  Us2=l  ’,  fOl 
WRITE  (10,*)  ’FT  at  HS  when  Usl=l  andUs2=0-:  ’,fl0 
WRITE  (10,*)  ’FT  at  HS  when  Usl=l  andUs2=I  ’,  fll 
WRITE  (10,*)  'Error  caused  by  Team  Su-aiegy  ’,  peteam 

DO50ic=  l.ib 
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WRITE  (1 1,1000)  TL(ic).  TU(ic).  (cbar(id,ic).id=l,10) 
50  CONTINUE 

DO70ie=  l.ib 

WRITE  (12,1000)  TL(ie),  TU(ie).  (cbar(ig.ie),ig=l  1,20) 
70  CONTINUE 

DO90ih=  1,21 
WRITE  (13,*)  ccc(ih),  cfmin(ih) 

90  CONTINUE 

DO  110  ii=  1,21 
WRITE  (14,*)  opthr(ii),  cfmin(ii) 
no  CONTINUE 

DO  130  ij  =  1,21 
WRITE  (15,*)  ccc(ij),  opthr(ij) 

130  CONTE'njE 

c  Format  statements 

1000  FORMAT  (’  ’T6.3,1X,H(F6.4,1X)) 

STOP 

ExND 


0%%%%  %%%%%%%%%%%%% 


SUBROUTINE  QfunE(xx,  erfcx) 


c  This  function  calculates  the  error  function  and  the  complimentary 
c  error  function  for  the  value  "xx" 


c  Accuracy  is  to  within  1.5E-07. 


REAL  X,  XX,  erfcx 

REAL  al,  a2,  a3,  a4,  a5,  p,  pi,  t 


pi  =  3.141592654 
X  =  ABS(xx) 


al  =  0.319381530 
a2  =  -0.356563782 
a3  =  1.781477937 


a4  =  -1.8..  1255978 
a5  =  1.330274429 

p  =  0.2316419 
t=  1.0/(1.0  +  p*x) 

si  =  al*i  +  a2’*i*’''2  +  a3*l**3  +  a4*t**4  -t-  a5*t**5 
s2  =  sl*EXP(-{x*»2)/2.0) 

IF  (XX  .GE.  0.0)  THEN 
erfcx  =  s2/SQRT(2.0*pi) 

ELSE  IF  (XX  LT.  0.0)  THEN 
erfcx  =  1.0  -  s2/SQRT(2.0*pi) 

ENT  IF 

RETLTLN 

END 


SUBROLTINE  FTNDMIN(array,isize,minindex) 

REAL  array(isize) 

INTEGER  minmdex 

uii  =  1 

1 1  CO^TL^^■E 

DO  10  1  =  int+1,  isize 
IF  (arrayOnt)  .GT.  arr3y(i))  THEN 
minmdex  =  i 
int  =  minindex 
GOTO  1 1 
END  IF 

10  CONTLN'L'E 


REILILN 

ENT) 


APPENDIX  C 


Program  for  Cost  Evaluation  of  2/3-Sensor-System 


The  program  listed  in  Appendix  C  is  used  /or  numerical  evaluation  of  2/3- 
Sensor-System’s  expected  costs  which  is  described  by  (4.2.1)  substituted  with  (4.5.1), 
(4.5.2),  (4.5.3),  (4.5.4),  and  (4.5.5). 

Figure  4.3,  Figure  4.4,  Figure  4.5,  Figure  4.6,  and  Table  5.3  are  constructed  by 
using  outputs  from  this  program. 
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PROGRAM  TW03SENSYS 


^••*»*** *•••«*«•*«•*•••**»•*«•**•••»*«••»**«»«*•*»•***•**«•*»•*«*•*****•*»*** 

c  Author.  Howard  C.  Choe 
c  Organization;  Department  of  Electrical  Engineering 
c  The  University  of  Virginia,  Charlottesville 

c  This  program  numerically  evaluates  e.xpected  costs  of 
c  a  ihree-scnsor-system  which  uses  team  strategies  for  Gaussian 
c  statistics. 

c  Variable  Description 

c  pO  ;  a  priori  probabilitv  of  "No  target  exists”,  HO 
c  pi  ;  a  priori  probabilii  /  of  "Targti  exists",  HI 
c  For  Host  Sensor 

c  mhO  :  mean  value  of  HO  received  by  HS 
c  mhl  ;  mean  value  of  HI  received  by  HS 
c  shO  :  standard  deviation  of  HO  received  by  HS 
c  shl  :  standard  deviation  of  HI  received  by  HS 
c  For  Slave  Sensor  i 

c  ms!0  :  mean  value  of  HO  received  by  SS  2 
c  msl!  ;  mean  value  of  HI  received  by  SSI 
c  sslO  :  standard  deviation  of  HO  received  by  SSI 
c  ssl  1  ;  standard  deviation  of  HI  received  by  SSI 

c  For  Slave  Sensor  2 

c  ms20  ;  mean  value  of  HO  received  by  SS2 

c  ms21  ;  mean  value  of  HI  received  by  SS2 

c  ss20  ;  standard  deviation  of  HO  received  by  SS2 

c  ss21  ;  standard  deviation  of  HI  received  by  SS2 

c  Thresholds 

c  For  Host  Sensor 

c  TL 1  ;  lower  tlircshold  1 

c  TL2  :  lower  threshold  2 

c  TUI  :  upper  threshold  1 

c  TU2  :  upper  threshold  2 

c  For  Slave  Sensor  1 

c  Tssl  ;  LRT  optimal  threshold 

c  For  Slave  Sensor  2 

c  Tss2  :  LRT  optimal  threshold 

c  Pre-cost  constant 

c  cOO  ;  deciding  HO  given  HO 

c  ciO  :  deciding  HI  given  HO 

c  cll  ;  deciding  Hi  gi-.cn  HI 

c  cOl  :  deciding  HO  given  HI 
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c  Team  effort  cost 

c  ccc  1  :  communication  cost  constant  for  communicating  wiih  one  sensor 

c  ccc2  :  communication  cost  constant  for  communicating  with  two  sensor 

c  Declaration 
REAL  pO,  pi 

REAL  mhO,  mhl,  mslO,  msl  1,  ms20,  ms21 
REAL  shO,  shl,  sslO,  ssll,  ss20,  ss2l 
REAL  cOO.clO,  clLcOl 
REAL  Tssl.  Tss2 

REAL  ccci(2l),  ccc2(21) 

REAL  cfmin(21),  opthr(21) 

REAL  TUl(lOOl).  TU2(l00l).  TLl(lOOl).  TL2(1001) 

REAL  pehs(lOOl),  pzcomUlOOl),  pzcom2(1001),  cfflOOl) 

REAL  cbar(2l. 1001) 


DATA  CCC2/0.00.  0.05,  0.10. 0.15. 0.20. 

•  0.25.0.30.0.35.0.40.0.45. 

•  0.50. 0.55.  0.60. 0.65. 0.70. 

•  0.75. 0.80. 0.85. 0.90. 0.95. 1.00/ 

c  Determination  of  communicauon  cost  constant  for  communicating 
c  wiih  2  slave  sensors 

D05  1=  1.21 
cccl(i)  =  0.5*ccc2(i) 

5  CONTINUE 

c  The  input  data 


p0  =  0.5 
pi  =0.5 

mhO  =  -1.0 
mhl  =  1.0 
mslO  =  -1.0 
msl  1  =  1.0 
ms20  =  -1.0 
ms21  =  1.0 

sigma  =  1.0 
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cOO  =  0.0 
clO=  1.0 
cl  1  =0.0 
cOl  =  1.0 

c  Assign  all  standard  deviation  to  the  same  value 

shO  =  sigma 
shl  =  sigma 
sslO  =  sigma 
ssll  =  sigma 
ss20  =  sigma 
ss21  =  sigma 

c  Evaluation  of  the  pre-caiculated  threshold  for  the  team  strategy  and 
c  for  the  slave  sensor 

plamssl  =■  (clO  -  c00)/(c01  -  cll) 
plamss2  =  (clO  -  c00)/(c0l  -  cl  1) 
plamt  =  plamssl 

c  Evaluation  of  the  ratio  between  a  priori  probabilities  and  LRT  threshold 
c  for  the  host  sensor  and  the  slave  sensor,  considering  each  sensor  is 
c  centralized  individually. 

tO  =  pO/pl 

Ths  =  (mhO  +  mhl)/2.0  +  (sigma**2/(mhl  -  mhO))*LOG(plamt*tO) 
Tssl  =  (mslO  +  msl  l)/2.0 

*  +  (sigma**2/(msll  -  mslO))*LOG(plamssl*iO) 

Tss2  =  (ms20  +  ms21)/2.0 

*  +  (sigma**2/(ms21  -  ms20))*LOG(plamss2*i0) 

c  Slave  sensor  1 

bslO  =  (Tssl  -  msl0)/ssl0 
bsll  =  (Tssl  -  msllVssll 

CALL  QfunE(bsl0,  QbslO) 

CALL  QfunE(bsll,Qbsll) 

c  Slave  sensor  2 

bs20  =  (Tss2  -  ms20)/ss20 
bs21  =  (Tss2  -  ms21)/ss21 
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CALL  QfunE(bs20,  Qbs20) 

CALL  QfunE(bs21,  Qbs21) 

c  FOR  COMMUNICATING  WITH  1  SLAVE  SENSOR 
c  Evaluation  of  the  final  threshold  after  communication 

fsO  =  plamssl*t0*(1.0-Qbsl0)/(1.0-Qbsll) 
fsl  =  plamssl*tO*QbslO/Qbsll 

c  Evaluation  of  Q(y)-funciion  values  using  the  above  values 

fsOhO  =  (fsO  -  mhOVshO 
fslhO  =  (fsl  -  mhOVshO 

CALL  QfunEffsOhO,  QfsOhO) 

CALL  QfunEffslhO,  QfslhO) 

fsOhl  =(fsO-mhl)/shl 
fslhl  =  (fsl  -  mhiyshl 

CALL  QfunE(fsOhl,  QfsOhl) 

CALL  QfunE(fslhl.Qfslhl) 

c  Calculation  of  an  error  probability  by  the  team  strategy 
c  with  communication  with  1  slave  sensor 

peteaml  =  (Qfs0h0*(1.0-Qbsl0)  +  QfslhO’QbslO)  *  pO 
•  +  (( 1 .0-QfsOh I )•(  1 .0-Qbs  1 1 )  +  (1 .0-Qfs Ih  1  )*Qbs  11 )  *  pi 

c  FOR  COMMUNICATING  WITH  2  SLAVE  SENSORS 
c  Evaluation  of  the  final  threshold  after  communication 


roo  =  plamss  1  •tO’f  1 .0-Qbsl0)*(l  .0-Qbs20)/((l  .0-Qbsl  1)*(  1 .0-Qbs2 1 )) 
IDI  =  plamssl*t0*(1.0-Qbsl0)*Qbs20/((1.0-Qbsll)*Qbs:i) 
no  =  plamss2*i0*Qbsl0*(1.0-Qbs20)/(Qosll*(1.0-Qbs21)} 
fll  =  plamss2*tO*QbslO*Qbs20/(Qbsl  l*Qbs21) 

c  Evaluation  of  Q(y)-function  values  using  the  above  values 

IDOhO  =  (fOO  -  mhOVshO 
roihO  =  (roi  -mhOVshO 
nOhO  =  (flO-mhO)/shO 
fllh0  =  (fll  -mhOVshO 
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CALL  QfunE(roOhO.  QfOOhO) 
CALL  QfunE(fOlhO.  QfDlhO) 
CALL  QfunE(flOhO.  QnOhO) 
CALL  QfunE(fl  IhO,  Qfl  IhO) 


IDOhl  =  (fOO  -  mhl)/shl 
roihl  =  (roi-mhl)/shl 
flOhl  =  (flO-mhl)/shl 
fllhl  =  (fll  -mhi)/shl 


CALL  QfunE(fOOhl.  QfDOhl) 
CALL  QfunE(fDlhl.  QfDlhl) 
CALL  QfunE(flOhl.  QflOhl) 
CALL  QfunE(fl  Ih  1.  Qfl  Ih  1) 


c  Calculation  of  an  error  probability  by  the  team  strategy 
c  with  communicating  with  2  slave  sensors 

peteam2  =  (  Qfl30h0*(L0<}bsl0)*(l.0-Qbs20) 

•  +Qroih0*(1.0-Qbsl0)*Qbs20 

•  +Qfl0hO*QbslO*(1.0-Qbs2O) 

•  +  Qf  1 1  hO*Qbs  10*Qbs20 )  *  pO 

•  +(  (I.0-Qro0hl)*(1.0-Qbsll)*(1.0-Qbs21) 

•  +(1.0-Qroihl)*(1.0-Qbsll)*Qbs21 

•  +(1.0-Qfl0hl)*Qbsll*(1.0-Qbs21) 

•  +(1.0-Qfllhl)*Qbsll*Qbs2l  )*pl 

c  Calculation  of  the  host  sensor  error  probability,  that  of 
c  communication  frequency  probability  with  1  sensor,  and  that 
c  of  communication  frequency  probability  with  2  sensors-pehs, 
c  pzcoml,  and  pzcom2-respectively. 

c  pehs  is  TLl  and  TUI  dependent,  pzcoml  depends  on  TLl  and  TL2, 
c  and  TU2  and  TUI,  and  pzcom2  is  dependent  upon  TL2  and  TU2. 


tine  =  0.02 
DO  10ia=  1,21 
ib  =  0 

DO  30  thr  =  Ths,  Ths  +  4.0,  tine 
ib  =  ib  +  1 
TUl(ib)  =  thr 

TLl(ib)  =  Ths  -  (FLOAT(ib)  -  1.0)*tinc 

TU2(ib)  =  (TUl(ib)-Ths)/2.0 
TL2(ib)  =  TLl(ib)+{Ths-TLl(ib))/2.0 


96 


c  Evaluation  for  pehs  for  various  TL  and  TU 

tulhO  =  (TUl(ib)  -  mhO)/shO 
Ulhl  =  (TLl(ib)-mhl)/shl 

CALL  QfunE(tulhO,  QtulhO) 

CALL  QfunEfUlhl,  QUlhl) 

pehs(ib)  =  QtulhO’pO  +  (1.0  -  QUlhl)*pl 

c  Evaluation  for  pzcoml  for  various  TLl,  TL2,  TUI,  and  TU2 

UlhO  =  (TLl(ib)  -  mhOyshO 
tl2h0  =  (TL2(ib)  -  mh0)/sh0 
tulhO  =  (TUl(ib)  -  mhOyshO 
tu2h0  =  (TU2(ib)  -  mhOyshO 

CALL  QfunE(tllhO,  QUlhO) 

CALL  QfunE(tl2hO,  Qtl2h0) 

CALL  QfunE(tulhO,  QtulhO) 

CALL  QfunE(tu2hO.  Qtu2h0) 

tllhl  » (TLl(ib)  •  mhl)/shl 
tl2hl  =  CIL2(ib)  ’  mhl)/shl 
tulhl  =  (TUl(ib)  •  mhl)/shl 
tu2hl  =(TU2(ib)-mhl)/shl 

CALL  QfunE(tllhl,  QUlhl) 

CALL  QfunE(U2hl,  QU2hl) 

CALL  QfunE(tuIhl,  Qtulhl) 

CALL  QfunE(tu2hl,  Qtu2hl) 

pzcoml(ib)  =  ( QUlhO  -  QU2hO  +  Qtu2h0  -  QtulhO )  •  pO 
•  +(QUlhl -Qtl2hl+Qtu2hl -Qtulhl) ’pi 


c  Evaluation  for  pzcom2  for  various  TL2  and  TU2 

pzcom2(ib)  =  (Qtl2h0  -  Qtu2h0)*p0  +  (QU2hl  -  Qtu2hl)*pl 
c  Evaluation  for  the  expected  probability  of  error  of  the  system,  COST 
cbar(ia,ib)  =  pchs(ib) 

•  +  (peteaml-pchs(ib)+cccl(ia))*pzcoml(ib) 

•  -*•  (pcteam2-pchs(ib)-»-ccc2(ia))*pzcom2(ib) 

•  +  (pehs(ib)-peteaml-peteam2-cccl(ia)-ccc2(ia)) 
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'pzcoml(ib)*pzcom2(ib) 


cf(ib)  =  cbar(ia,ib) 

30  CONTINUE 

c  Extract  the  minimum  system  cost  for  each  case  of  ccc(ia) 

CALL  FINDMIN(cfab4nini) 
cfmin(ia)  =  cf(mini) 
opthr(ia)  =  TUl(mim) 

10  CONTINUE 

c  Write  OUTPUT  DATA . 

WRITE  (10,*)  ’Ratio  of  a  priori  probabilities  lO 


WRITE  (10.*)  ’LRT  Threshold  of  SSI - Tssl 

WRITE  (10,*)  *LRT  Threshold  of  SS2 - Tss2 

WRITE  (10,*)  ’Eaor  caused  by  TS  using  SSI  only peteaml 

WRITE  (10,*)  'FT  at  HS  when  Usl=0 - fsO 

WRITE  (10  •)  ’FT  at  HS  when  Usl=l - fsl 

WRITE  (10.*)  ’  ’ 

WRITE  (10.*)  ’Enor  caused  by  TS  using  SSI  &  SS2:’,  peteam2 


WRITE  (10.*)  ’FT  at  HS  when  Usl=0  and  Us2=0  — fOO 
WRITE  (10,*)  ’FT  at  HS  when  Usl=0  and  Us2=l  — fOl 
WRITE  (10,*)  ’FT  at  HS  when  Usl=l  and  Us2=0  — flO 
WRITE  (10,*)  ’FT  at  HS  when  Usl=l  and  Us2=l  — fl  1 

D045ja=  Lib 

WRITE  (9.*)  TLl(ja).  TL2(ja),  TU2(ja),  TUlOa) 

45  CONTINUE 

DO  50  ic=  Lib 

WRITE  (11.1000)  TUI(ic).  (cbar(id.ic).id=L10) 

50  CONTINUE 

DO70ie»  Lib 

WRITE  (12.1000)  TUl(ic),  (cbar(ig.ie).ig=l  1,20) 

70  CONTINUE 

DO90ih=  1,21 

WRITE  (13.*)  ccc2(ihi,  cfmin(ih) 

90  CONTINLT 
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DO  110u=  1,21 
WRITE  (14,*)  opthr(ii),  cfmin(ii) 
110  CONTINUE 

DO  130  ij  =  1,21 
WRITE  (15.*)  ccc2(ij)  opihr(ij) 
130  CONTINUE 

c  Format  statements 

1000  FORMAT(’ M!(F6.4.1X)) 

STOP 

END 


SUBROUTINE  QfunE(xx,  erfcx) 

c  This  function  calculates  the  error  function  and  the  complimentary 
c  error  function  for  the  value  "xx" 

c  Accuracy  is  to  within  1.5E-07. 

REAL  X,  XX,  erfcx 

REAL  al,  a2,  a3,  a4,  a5,  p,  pi,  t 

pi  =  3.141592654 
X  =  ABS(xx) 

al  =  0.319381530 
a2  =  -0.356563782 
a3  =  1.781477937 
a4  =  -1.821255978 
a5  =  1.330274429 

p  =  0.2316419 
t=  1.0/(1.0-t-p*x) 

si  =  al*t  +  a2*t**2  +  a3*t**3  +  a4*t**4  -t-  a5*t**5 
s2  =  sl*EXP(-(x**2)/2.0) 

IF  (XX  .GE.  0.0)  THEN 
erfcx  =  s2/SQRT(2.0*pi) 


99 


ELSE  IF  (XX  J-T.  0.0)  THEN 
erfcx  =  1.0  -  s2/SQRT(2.0*pi) 
END  IF 

RETURN 

END 


SUBROUTINE  FINDMIN(array,isize,minindex) 

real  array(isize) 

INTEGER  minindex 


inta  1 

11  CONTINUE 

DO  10i  =  int+l,isize 
IF  (array(im)  .GT.  array(i))  THEN 
minindex  =  i 
int  =  minindex 
GOTO  11 
END  IF 

10  CONTINUE 

RETURN 

END 
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APPENDIX  D 


Program  for  Calculation  of  Dubious  Decision  Probabilities 


In  this  appendix,  a  program  UNPRO  (UNcenainty  PRObability)  is  attached. 
This  program  evaluates  the  dubious  decision  probability  at  the  host  sensor  in  2SS, 
3SS,  and  23SS  when  the  optimal  thresholds  for  given  communication  cost  constants 
are  known.  These  thresholds  can  be  obtained  from  the  programs  attached  in  Appen¬ 
dix  A,  Appendix  B,  and  Appendix  C. 

The  information  obtained  by  this  program  are  plotted  in  Figure  5.1,  Figure  5.2, 
and  Figure  5.3.  The  tabulated  data  can  also  be  found  in  Table  5.8,  Table  5.9,  and 
Table  5.10. 
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PROGRAM  UNPRO 


^* •«•»•*••* *««»»«*•*•«******«»****%*«•»«*••****•*»»»«»•*•«•*»» •»**•»»•»***•» 

c  This  program  UNPRO  (UNcenainty  PRObability)  is  written  to  evaluate 
c  the  probability  of  the  observation  that  falls  in  the  uncertainty 
c  region  of  the  host  sensor.  This  program  reads  in  the  optimal  threshold 
c  locations  which  are  evaluated  using  programs  such  as  2sensys.f, 
c  Ssensys.f,  and  2/3sensys.f  (These  programs  are  listed  in  Appendix  A,  B, 
c  and  C,  respectively.). 

^•«««»*****»«***«*»««*«*«i|i***««»«««»»«»«»«*»»»*«»»***«*4i»«*»*»*»»*»**««»**«» 

c  Declaration 

REAL  otp(3.21).error(3.21) 

REAL  qlimR(3.21),  qlimL(3^1) 

REAL  ytintR(3,21),  ytintL(3^1) 

REAL  prob(3^1) 

c  Open  data  file  and  rewind 

OPEN  (UNIT=llJTLE=fortir) 

OPEN  (UNIT=12JTLE=’foal2’) 

OPEN  (UNTT=13.FILE=’fortl3’) 

REWIND  (11) 

REWIND  (12) 

REWIND  (13) 

c  Statistics  for  Gaussian  observation 
avg  =  1.0 

sigma  =  1.0 

c  Read  in  input  data  from  files  opened 
DO  10i=  1,21 

READ  (11,*)  otp(l,i),  error(l,i) 

READ  (12,*)  otp(2,i),  error(2,i) 

READ  (13  •)  otp(3.i).  error(3,i) 

10  CONTINUE 

c  Calculation  of  limits  (qlimR  and  qlimL)  for  Ofy)-funciion 
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c  and  evaluate  the  corresponding  probabilities. 


DO  30  i*  1.3 
DO  50  j=  1,21 

qlimR(i  j)  =  (otp(i,j)-avg)/sigma 
qlimL(i.j)  =  (-otp(iJ)-avg)/sigma 

xxR  =  qliiiiR(i,j) 

CALL  QfunE(xxR,erfcxR) 

xxL  =  qlimLCi.ji 
CALL  QfunE(xxL,erfcxL) 

ylintR(i  j)  =  erfcxR 
ytintL(i,i)  =  erfcxL 
50  CONTINUE 
30  CONTINUE 

c  Evaluate  the  probability  of  dubious  decision  for  one  hypothesis 

DO  70  i*  1,3 
DO90j=  1,21 

prob(ij)  =  ytintL(ij)  -  ytiniR(ij) 

90  CONTINUE 
70  CONTINUE 

c  prob(,)  is  multiplied  by  2  J  since  prob(,)  is  the  probability  of 
c  a  observation  that  la.^ds  in  the  uncertainty  region  of  HS  given  HI. 
c  The  dubious  decision  probability,  when  HO  is  considered,  is  the  same. 

DO  110  i^  1,21 

WRITE  (21,1000)  otp(l,i).  prob(l  i)*2.0 
WRITE  (22,1000)  otp(2,i).  prob(2,i)*2.0 
WRITE  (23,1000)  oip(3.i).  prob(3,i)*2.0 
110  CONTINU'E 

STOP 

1000  FORMAT  ('  ’,F8.3,1X.F8.4) 


END 
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SUBROUTINE  QfunE(xx,  erfcx) 


c  This  function  calculates  the  error  function  and  the  complimentary 
c  error  function  for  the  value  "xx" 

c  Accuracy  is  to  within  1.5E-07. 

REAL  X,  XX,  erfcx 

REAL  al,  a2,  a3,  a4,  a5,  p,  pi.  t 

pi  =  3.141592654 
X  =  ABS(xx) 

al  =  0.319381530 
a2  =  -0.356563782 
a3  =  1.781477937 
a4  =  -1.821255978 
a5  =  1.330274429 

p  =  0.2316419 
t=  1.0/(1.0  +  p*x) 

si  =  al’t  +  a2*t**2  +  a3*t**3  +  a4»t*»4  +  a5*t‘'*5 
s2  =  sl*EXP(-(x**2)/2.0) 

IF  (XX  .GE.  0.0)  THEN 
erfcx  =  s2/SQRT(2.0*pi) 

ELSE  IF  (XX  .LT.  0.0)  THEN 
erfcx  =  1.0  -  s2/SQRT(2.0*pi) 

END  IF 

RETURN 

END 


APPENDIX  E 


Program  Listing  of  System  Simulation 


The  program,  SENSIM  (SENsor  SIMuIation),  simulates  2SS,  3SS,  and  2/3SS 
for  Gaussian  observations.  This  program  incorporates  the  programs  listed  in  Appen¬ 
dix  A,  B,  and  C.  These  programs  -  TWOSENSYM,  THREESENSYM,  and 
TW03SENSYM  -  become  subroutines  named  SETBAND2,  SETBAND3,  and  SET- 
BAND23,  respectively.  These  subroutines  return  the  optimal  thresholds  location  and 
the  final  decision  thresholds  at  the  host  sensor  for  a  given  communication  cost  con¬ 
stant.  There  are  other  subroutines  which  are  used  to  generate  Gaussian  random 
observation  (or  Gaussian  random  number).  Because  of  the  mutually  independent 
observations  among  sensors,  each  sensor  is  provided  with  its  own  Gaussian  random 
observation  generator. 

Outputs  from  this  simulation  are  presented  in  Figure  6.1,  Figure  6.2,  and  Figure 
6.3.  Tabulated  data  of  these  figures  are  in  Table  6.1  and  Table  6.2. 

In  the  subroutine  SETBAND2,  SETBAND3,  and  SETBAND23,  all  the  com¬ 
ments  are  omitted  since  they  are  the  same  as  the  programs  attached  in  Appendix  A, 
B,  and  C. 
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PROGRAM  SENS  EM 


c  AUTHOR:  HOWARD  C.  CHOE 
c  ORGANIZATION:  Department  of  Elcrirical  Engineering 
c  University  of  Virginia,  CharlottcsvilJfr 

c  This  program  simulates  the  sensor  systems  ( 2. 3.  and  2/3  sensor 
c  system)  which  uses  team  strategies.  The  host  sensor  and  the  slave 
c  sensors  receive  independent  observations  from  the  binary  hypothesis 
c  environment  under  Gaussian  model. 

c  Find  the  optimum  thresholds  of  the  host  sensor  (TL  &  TU  or  TLl,  TL2, 
c  TU2,  and  TUI)  for  the  different  system. 

VvT?ITE  (6,*)  ’ENTER  ccc  for  each  system,  2, 3.  &  23’ 

READ  (5,*)  ccc2,  ccc3,  ccc23 

WRITE  (6,*)  ’ENTER  #  of  iterations  desired’ 

READ  (5,*)  nter 

WRITE  (6.*)  ’ENTER  #  seed  for  a  random  #  generation’ 

READ  (5,*)  iseed 


TLRT  =  0.0 

CALL  SETBAND2(ccc2,TL2,TU2E20T21,T2) 

CALL  SETB AND3(ccc3 .TL3,TU3E300E301  J310,n  1 1  ,T3 1  ,T32) 
CALL  SETB AND23(ccc23,TL3 1  ,TL32,TU32,TU3 1 , 

■'  F0,F1E00E01E10.F11,T231,T232) 


WRITE  (10,*)  '*••••••••  ONE-SENSOP -SYSTEM 

WRITE  (10,*)’  ’ 

WRITE  (10,*)  ’LRT  Threshold  — . :  ’,TLRT 

WRITE  (10,*)  ’•••••••••  TWO-SENSOR-SYSTEM  •••••••••’ 

WRITE  (10,*)  ’Lower  Threshold - ;  TL2 

WRITE  (10,*)  ’Upper  Threshold - :  ’,1X12 


WRITE  (10,*)  ’LRT  threshold  of  Slave  Sensor  ’,  T2 
WRITE  (10,*)  ’Final  Threshold  when  Us  =  0  :  ’,  F20 
WRITE  (10,*)  ’Final  Threshold  when  Us  =  I  :  ’,  F21 
WRITE  (10.*)’  ’ 

WRITE  (10,*)  ■**••••>■•  THREE-SENSOR-SYSTEM********’ 
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WRITE  (10, •)  ’Lower  Threshold - :  ’.TLS 

WRITE  (10,»)  ’Upper  Threshold - :  ’.  TU3 


WRITE  (10,*)  ’LRT  threshold  of  SSI - :  ’.T3I 

WRITE  (10.*)  ’LRT  threshold  of  SS2 - :  T32 

WRITE  (10.*)  ’FT  when  Usl  =  0  &  Us2  =  0  F300 

WRITE  (10,*)  ’FT  when  Usl  =  0  &  Us2  =  1  — :  F301 

WRITE  (10.*)  ’FT  when  Usl  =  1  &Us2  =  0--:  ’.F310 
WRITE  (10,*)  ’FT  when  Usl  =  I  &  Us2  =  1  — :  ’.  F311 
WRITE  (10.*)’  ’ 

WRITE  (10.*)  *******  TWO/THREE-SENSOR-SYSTEM  ******’ 


WRITE  (10.*)  ’Lower  Threshold  1 - :  TL31 

WRITE  (10,*)  ’Lower  Threshold  2 - ;  ’,  TL32 

WRITE  (10.*)  ’Upper  Threshold  2 - :  ’.  TU32 

WRITE  (10,*)  ’Upper  Threshold  1 - :  ’,TU31 

WRITE  (10.*)  ’LRT  threshold  of  SSI - :  ’,  T231 

WRITE  (10.*)  ’LRT  threshold  of  SS2 - :  ’.  T232 

WRITE  (10.*)  ’FT  when  Usl  =  0 - :  ’.  FO 

WRITE  (10.*)  ’FT  when  Usl  =  0 - : FI 


WRITE  (10.*)  ’FT  when  Usl  =  0  &  Us2  =  0  FOO 
WRITE  (10.*)  ’FT  when  Usl  =  0  &  Us2  =  1  ’.  FOl 

WRITE  (10,*)  ’FT  when  Usl  =  1  &  Us2  *  0  — :  ’.  FIO 
WRITE  (10.*)  ’FT  when  Usl  =  1  &  Us2  =  1  — :  ’.  FI  1 

c  Standard  deviation  of  observation 

sigma  =  1.0 

c  ITERATION  STARTS 


ieO  =  0 
iel  =  0 

icdl  =  0 
ifal  =  0 
imtl  =  0 

icd2  =  0 
ifa2^0 
imt2  =  0 


icd3  =  0 
ifa3  =  0 
imi3  =  0 

icd23  =  0 
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ifa23  =  0 
imt23  a  0 

c  Get  system  clock  time  for  random  seeds 

c  itm  =  mciockO 

c  WRITE  (6,*)  ’itm  =  itm 

c  iseed  =  74591  +  2*MOD(1000*itm.500) 

CALL  SRAND(isced) 

DO  10  ia  =  1,  ntcr 

c  Generate  Environment 

11  CALL  GENENV(ienv) 
env  =  FLOAT(ienv) 

IF  (env  EQ.  -1.0)  THEN 
ieO  =  ieO  +  1 

ELSE  IF  (env  EQ.  1.0)  THEN 
iel  =  iel  +  1 
ELSE 

WRITE  (6  •)  ’###  Generated  ENV  is  NOT  either  -1  or  1  ###’ 
GOTO  11 
END  IF 

c  Generate  Observations  at  each  sensors 
c  For  1-Sensor-System 

CALL  HSl(env.sigma.yhl) 

c  For  2-Sensor-System 
c  CALL  HS2(env,sigma,yh2) 

CALL  SS2(env4igma.ys2) 

c  For  3-Sensor-Systcm 
c  CALL  HS3(env,sigma,yh3) 

c  CALL  SS31(env, sigma, ys3l) 

CALL  SS32(cnv,sigma,ys32) 

c  For  23-Scnsor-Systcm 
c  CALL  HS23(cnv,sigma,yh23) 

c  CALL  SS231(cnv,sigma,ys231) 

c  CALL  SS232(env,sigma,ys232) 

yh2  =yhl 
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yh3  =yhl 
yh23  =yhl 

ys3 1  =  ys2 
ys231  =  ys2 
ys232  =  ys32 

c  WRITE  (33,1 100)  env,yhl,yh2,ys2,yh3,ys31,ys32,yh23,ys231,ys232 
cllOO  FORMAT  (”.10(F7.3, IX)) 

c  Single  Sensor  using  LRT  Threshold 
IF  (yhl  iE.  TLRT)  THEN 
Uhl  =  -1.0 

ELSE  IF  (yhl  .GT.  TLRT)  THEN 
Uhl  =  1.0 
END  IF 

iuhl  =  INT(uhl) 

c  Count  False  alarm,  missing  target,  and  correct  detection 
IF  (iuhl  .EQ.  ienv)  THEN 
icdl  =  icdl  +  1 

ELSE  IF  (iuhl,EQ.l  .AND.  ienv£Q.-l)  THEN 
ifal  =  ifal  +  1 

ELSE  IF  (iuhl£Q.-l  .AND.  ienv.EQ.l)  THEN 
imtl  =  imtl  +  1 
END  IF 

c  Use  Team  Strategies  to  make  Decision 
c  For  2-Sensor-Sysiem 

IF  (yh2  .LE.  'iL2)  THEN 
uh2  =-1.0 

ELSE  IF  (yh2  .GE.  TU2)  THEN 
uh2  =  1.0 

ELSE  IF  (yh2.GT.TL2  .AND.  yh2.LT.TU2)  THEN 
IF  (ys2  .LE.  T2)  THEN 
us2  =  -1.0 

IF  (yh2  .LE.  F20)  THEN 
uh2  =  -1.0 

ELSE  IF  (yh2  .GT.  F20)  THEN 
uh2  =  1.0 
END  IF 

ELSE  IF  (ys2  .GT.  T2)  THEN 
"c2  =  1.0 

IF(yh2.LE.  F21)THEN 
uh2  =  -1.0 
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ELSE  IF  (yh2  .GT.  F21)  THEN 
uh2  =  1.0 
END  IF 
END  IF 
END  IF 

iuh2  =  INT(uh2) 

c  Count  False  alann,  missing  target,  and  correct  detection 
IF  (iuh2  .EQ.  ienv)  THEN 
icd2  =  icd2  +  1 

ELSE  IF  (iuh2£Q.l  .AND.  ienvEQ.-l)  THEN 
ifa2  =  ifa2  +  1 

ELSE  IF  (iuh2EQ..l  .AND.  ienv.EQ.l)  THEN 
imt2  =  imt2  +  1 
END  IF 

c  For  3-Se''sor-System 
rr(yh3  EE.  TL3)THEN 
uh3  =-1.0 

ELSE  IF  (yh3  .GE.  TU3)  THEN 
uh3  =  1.0 

EI,.SE  IF  (yh3.GT.TL3  .AND  yh3.LT.TU3)  THEN 
IF  (ys3l.LE.T31  .AND.  ys32.LE.T32)  THEN 
us31  =  -1.0 
us32  =  -1.0 

IF  (yh3  .LE.  F300)  THEN 
uh3  =  -1.0 

ELSE  IF  (yh3  .GT.  F300)  THEN 
uh3  =  1.0 
END  IF 

ELSEIF(ys31EE.T3l  .AND.  ys32.GT.T32)  THEN 
us31  =  -1.0 
us32=  1.0 

IF  (yh3  .LE.  HOI)  THEN 
uh3  =  -1.0 

ELSE  IF(yh3  .GT.  F301)  THEN 
uh3  =  ’  0 
END  IF 

ELSE  IF  (ys31.GT.T31  .AND.  ys32.LE.T32)  THEN 
us31  =  1.0 
us32  =  -1.0 

IF(yh3  .LE.F310)  THEN 
uh3  =  -1.0 

ELSE  IF  (yh3  .GT.  F310)  THEN 
uh3  =  1.0 
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END  IF 

ELSE  IF  (ys31.GT.T31  .AND.  ys32.GT.T32)  THEN 
us3l  =  1.0 
us32=  1.0 

IF(yh3  .LE.F311)THEN 
uh3  =  -1.0 

ELSE  IF  (yh3  .GT.  F3 1 1)  THEN 
uh3  =  1.0 
END  IF 
END  IF 
END  IF 

iuh3  =  INT(uh3) 

c  Count  False  alarm,  missing  target,  and  correct  detection 
IF  (iuh3  .EQ.  ienv)  THEN 
icd3  =  icd3  +  1 

ELSE  IF  (iuh3.EQ.l  .AND.  ienv.EQ.-l)  THEN 
ifa3  =  ifa3  +  1 

ELSE  IF  (iuh3£Q.-l  .AND.  ienv.EQ.l)  THEN 
imt3  =  imt3  +  1 
END  IF 

c  For  2/3-Sensor-System 

IF(yh23  IE.  TL31)THEN 
uh23  =  -1.0 

ELSE  IF  (yh23  .GE.  TU31)  THEN 
uh23  =  1.0 

ELSE  IF  (yh23.GT.TL31  .AND.  yh23.LT.TL32  .OR. 

•  yh23.GT.TU32  .AND.  yh23.LT.TU3 1)  THEN 
IF(ys231  .LE.  T231)THEN 
US231  =  -1.0 
IF  (yh23  .LE.  FO)  THEN 
uh23  =  -1.0 

ELSE  IF  (yh23  .GT.  FO)  THEN 
uh23  =  1.0 
END  rr 

ELSE  IF  (ys23  i  .GT.  T231)  THEN 
us231  =  1.0 

IF(yh23  .LE.  FI)  THEN 
uh23  =  -1.0 

ELSE  IF  (yh23  .GT.  FI)  THEN 
uh23  =  1.0 
END  IF 
END  IF 

ELSE  IF  (yh23.GE.TL32  .AND.  yh23.LE.TU32)  THEN 
IF  (ys231.LE.T231  .AND.  ys232.LE.T232)  THEN 
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us23l  =  -1.0 
us232  =  -1.0 

IF  (yh23  .LE.  FOO)  THEN 
uh23  =  -1.0 

ELSE  IF  Ot.23  .GT.  FOO)  THEM 
uh23  =  1.0 
END  IF 

ELSE  IF(ys231iE.T231  ..4ND.  ys232.GT.T232)  THEN 
us231  =  -1.0 
us232=  1.0 

IF  (yh23  .LE.  FOl)  THEN 
uh23  =  -1.0 

ELSE  IF  (yh23  .GT.  FOl)  THEN 
uh23  =  1.0 
END  IF 

ELSE  IF(ys23I.GT.T23I  .AND.  ys232.LE.T232)  THEN 
us231  =  1.0 
us232  =  -1.0 

IF  (yh23  .LE.  FIO)  THEN 
uh23  =  -1.0 

ELSE  IF  (yh23  .GT.  FIO)  THEN 
uh23  »  l.O 
ENDEF 

ELSE  IF  (ys231.GT.T231  .AND.  ys232.GT.T232)  THEN 
us231  =  1.0 
us232=  1.0 

IF  (yh23  .LE.FIDTHEN 
uh23  =  -1.0 

ELSE  IF  (yh23  .GT.  FI  1)  THEN 
uh23  =  1.0 
END  IF 
END  IF 
END  IF 

iuh23  =  INT(uh23) 

c  Count  False  alann,  missing  target,  and  correct  detection 
IF  (iuh23  .EQ.  ienv)  THEN 
icd23  s  icd23  +  1 

ELSE  IF  (iuh23.EQ.l  .AND.  ienv.EQ.-l)  THEN 
ifa23  =  ifa23  +  1 

ELSE  IF  (iuh23.EQ.-l  .AND.  ienv.EQ.l)  THEN 
iml23  =  imt23  +  1 
END  IF 

10  CONTINUE 
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c  Find  the  percentage  of  O’s  and  I’s  in  total  environment  generated 


pevO  =  100.0*FLOAT(ie0)/FLOAT(ia) 
pevl  =  100.0*FLOAT(iel)/FLOAT(ia) 

WRITE  (lO.*)  ’  ■ 

WRITE  (10,*)  ’@@@@@@  %  of  0  or  1  of  the  Environment 
WRITE  (10.*) '%  of  Os  pevO 
WRITE  (10.*)  ’%  of  Is  pevl 

c  Find  the  percentage  of  correct  detection,  faise  alarm,  and  missing  target 
c  for  each  system. 

c  1 -Sensor-System 

pcdl  =  100.0*FLOAT(icdl)/FLOAT(ia) 
pfal  =  100.0*FLOAT(ifal)/FLOAT(ia) 
pmtl  =  IOO.O*FLOAT(imtI)/FLOAT(ia) 

c  2-Sensor-System 

pcd2  =  100.0*FLOAT(icd2)/FLOAT(ia) 
pfa2  =  100.0*FLOAT(ifa2)/FLOAT(ia) 
pmt2  =  100.0*FLOAT(imt2)/FLOAT(ia) 

c  3-Sensor-Sysiem 

pcd3  =  100.0*FLOAT(icd3)/FLOAT(ia) 
pfa3  =  100.0*FLOAT(ifa3)/FLOAT(ia) 
pmG  =  100.0*  FLO  AT(imt3)/FLOAT(ia) 

c  23-Sensor-System 

pcd23  =  10C.O*FLOAT(icd23)/FLOAT(ia) 
pfa23  =  100.0*FLOAT(ifa23)/FLOAT(ia) 
pmt23  =  100.0*FLOAT(imt23)/FLOAT(ia) 

c  WRITE  the  percentages 
WRITE  (10,*)’  ’ 

WRITE  (10,*)  ’//  %  of  CD,  FA,  and  MT  for  1 -Sensor-System  V 

WRITE  (10,*)  ’Correct  Decision  %:  ’,  pcdl 

WRITE  (10,*)  ’False  Alarm  %  — ;  ’,  pfal 

WRITE  (10,*)  ’Missing  Target  %  :  ’,  pmtl 

WRITE  (10,*)  ’CD  -I-  FA  -t-  MT  in  %  :  pcdl-fpfal-t-pmtl 

WRITE  (10,*)’  ’ 

WRITE  (10,*)  ’//  %  of  CD,  FA,  and  MT  for  2-Scnsor-System\’ 
'VRITE  (10,*)  ’Correct  Decision  %:  ’,  pcd2 
WRITE  (10,*)  ’False  Alarm  %  - — :  ’,  pfa2 
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WRITE  (10,*)  ’Missing  Target  %  :  ptnt2 

WRITE  (10.*)  ’CD  +  FA  +  MT  in  %  :  pcd2+pfa2+pmt2 

WRITE  (10.*)’  ’ 

WRITE  (10,*)  ’//  %  of  CD,  FA,  and  MT  for  3-Sensor-System  V 

WRITE  (10,*)  ’Correct  Decision  %;  ’,  pcd3 

WRITE  (10.*)  ’False  Alarm  %  — :  ’,  pfa3 

WRITE  (10,*)  ’Missing  Target  %  :  ’,  pmt3 

WRITE  (10,*)  ’CD  +  FA  +  MT  in  %  :  ’,  pcd3+pfa3+pmt3 

WRITE  (10,*)’  ’ 

WRITE  (10,*)  ’//  %  of  CD,  FA,  and  MT  for  23-Sensor-System  V 

WRITE  (10,*)  ’Correct  Decision  %:  ’,  pcd23 

WRITE  (10.*)  ’False  Alarm  %  — :  pfa23 

WRITE  (10,*)  ’Missing  Target  %  :  ’,  pmt23 

WRITE  (10.*)  ’CD  +  FA  +  MT  in  %  :  ’.  pcd23+pfa23+pmt23 

STOP 

END 

Subroutine  for  the  generation  of  environment 
SUBROUTINE  GENENV(ienv) 


im  =  irandO 

xm  =  FLOAT(ifn)/FLOAT(2**15  -  1) 
renv  =  xm  -  0.5 

IF  (renv  IE.  0.0)  THEN 
ienv  =  -1 

ELSE  IF  (renv  .GT.  0.0)  THEN 
ienv  s  1 
END  IF 

RETURN 

END 


e'7o‘?o  Subroutine  for  Singlc-Scnsor-Systcm 


'7c‘7g9o?o%%%%%%%%%%<7o%  %  ?c  %  %  ?o  <7c  %  %  ?c  %  "c  '7c 


SUBROUTINE  HSl(env,sigma,yhl) 


a  =  0.0 

DO  10  i=  1,  12 
im  =  irandO 

xm  =  FLOAT(itn)/FLOAT(2**I5  -  I) 
a  =  a  +  xm 
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10  CONTINUE 


yhi  =  (a  -  6.0)*signia  +  env 

RETURN 

END 

c%%  Subroutines  for  2-Sensor-Systeni 
SUBROUTINE  HS2(env,sigma,yh2) 
a  =  0.0 

DO  10i=  1, 12 
im  =  irandO 

xm  =  FLOAT(im)/FLOAT(2**15  -  1) 
a  =  a  +  xm 
10  CONTINUE 

yh2  =  (a  -  6.0)*sigma  +  cnv 

RETURN 

END 


SUBROUTINE  SS2(env,sigma.ys2) 
a  =  0.0 

DO  10i=  1.  12 
im  =  irandO 

xm  =  FLOAT(im)/FLOAT(2**15  - 1) 
a  =  a  +  xm 
10  CONTINUE 

ys2  =  (a  -  6.0)*sigma  +  env 

RETURN 

END 

c%%  Subroutines  for  3-Scnsor-Systcm 
SUBROUTINE  KS3(env,sigma,yh3) 
a  =  0.0 

DO  10i=  1, 12 
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im  =  irandO 

xm  =  FLOAT(ini)/FLOAT(2**15  - 1) 
a  =  a  +  xm 
10  CONTINUE 

yh3  =  (a  -  6.0)*sigma  +  env 

RETURN 

END 


SUBROUTINE  SS31(env4igma,ys3I) 
a  =  0.0 

DO  10  i  =  1.  12 
im  =  irandO 

xm  =  FLOAT(im)/FLOAT(2**15  - 1) 
a=a+xm 
10  CONTINUE 

ys31  =  (a  •  6.0)*sigma  +  env 

RETURN 

END 


SUBROUTINE  SS32(env,sigma,ys32) 

a  =  0.0  • 

DO  10  i=  1,12 
im  =  irandO 

xm  =  FLOAT(im)/FLOAT(2**15  - 1) 
a  =  a  +  xm 
10  CONTINUE 

ys32  =  (a  -  6.0)*signia  +  env 

RETURN 

END 

c^o%  Subroutines  for  23-Scnsor-Sysicm  %7o‘7c%'7o9o%9o%%9o%%%%'7o%%^c'^c^c%'7c%?c?i 
SUBROUTLNE  HS23(cnv,sigma,yh23) 
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a  =  0.0 

DO  10  i=  1, 12 
im  =  irandO 

xm  =  FLOAT(im)/FLOAT('2**15  -  1) 
a  =  a  +  xm 
10  CONTINUE 

yh23  =  (a  -  6.0)*sigma  +  env 

RETURN 

END 


c 


SUBROUTINE  SS231(env,sigma,ys231) 


a  =  0.0 

DO  10i=  1,  12 
im  =  ixandO 

xm  =  FLOAT(im)/FLOAT(2**I5  -  1) 
a=a+xm 
10  CONTINUE 

ys23 1  =  (a  -  6.0)*sigma  +  env 

RETURN 

END 


c 


SUBROUTINE  SS232(cnv,sigma,ys232) 


a  =  0.0 

DO  10i=  1.  12 
im  =  irandO 

xm  =  FLOAT(im)/FLOAT(2**15  -  1) 
a  =  a  +  xm 
10  CONTINUE 

ys232  =  (a  -  6.0)*sigma  +  env 

RETURN 

END 
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SUBROUTINE  SETBAND2(ccc2.TL2.TU2.F20.F2I.T2) 
REAL  pO,  pi 

REAL  mhO,  mhl,  msO,  msl 
REAL  shO,  shl,  ssO,  ssl 
REALcOO.clO.cll.cOl 
REALTss 

REAL  ccc2,  cfmin,  oplthr,  opulhr 
REAL  TU(lOOl).  TL(lOOl) 

REAL  pehs(lOOl),  pzcom(lOOl) 

REAL  cbar(100l) 


p0  =  0.5 
pi  =0.5 


mhO  =  -1.0 
mhl  =  1.0 
msO  =  -l.O 
msl  =  1,0 

sigma  =  1.0 

cOO  =  0.0 
clO=  1.0 
cl  1  =0,0 
cOl  =  1.0 

shO  =  sigma 
shl  =  sigma 
ssO  =  sigma 
ssl  =  sigma 

plamss  =  (clO  -  c00)/(c01  -  cll) 
plamt  =  plamss 

lO  =pO/pl 

Ths  =  (mhO  +  mhl)/2.0  +  (sigma**2/(mhl  -  mhO))*LOG(plamt*iO) 
Tss  =  (msO  +  msl)/2.0  +  (sigma* *2/(msl  -  msO))*LOG(plamss*tO) 

aOs  =  (Tss  -  msO)/ssO 
als  =  (Tss  -  msl)/ssl 
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CALL  QfunE(aOs.  QaOs) 

CALL  QfunE(als,  Qals) 

fusO  =  plamss  •  lO  *  (1.0  -  Qa0s)/(1.0  -  Qals) 
fusl  =  plamss  *  lO  *  QaOs/Qals 

fusOhO  =  (fusO  -  mhO)/shO 
fuslhO  =  (fusl  -  mhO)/shO 
fusOhl  =  (fusO  -  mhl)/shl 
fuslhl  =  (fusl  -  mhl)/shl 

CALL  QfunE(fusOhO,  QfusOhO) 

CALL  QfunE(fuslhO.  QfuslhO) 

CALL  QfunE(fusOhl.  QfusOhl) 

CALL  QfunE(fuslhl,  Qfuslhl) 

peteam  =  (QfusOhO*(LO-QaOs)  +  QfuslhO’QaOs)  *  pO 
•  +  ((1.0-Qfus0hl)*(1.0-Qals)  +  (1.0-Qfuslhi)*Qals)  *  pi 

tine  =  0.02 
ib  =  0 

IX)  30  ihr  =  Ths.  Ths  +  2.0,  tine 
ib  =  ib  +  1 
TL'(ib)  =  ihr 

TL(ib)  =  Ths  -  (FLOAT(ib)  •  1.0)*iinc 

tuhO  =  (TL'(ib)  -  mhO)/shO 
ilhl  =  (TL(ib)  -  mhl)/shl 

CALL  QfunE(tuhO,  QtuhO) 

CALL  QfunEfUhl.QUhl) 

pchs(ib)  =  QtuhO*pO  +  (1.0  -  QiJhl)*pl 

tlhO  =  (TL(ib)  -  mhO)/shO 
tuhl  =(TU(ib)-mhl)/shI 

CALL  QfunE(tihO,  QtihO) 

CALL  QfunEduhl.Qluh!) 

pzcom(ib)  =  (Qtlhn  -  CtuhO)*pO  +  (Qtlhl  -  Qiuhl)»pl 

cbai-(ib)  =  pehs(ib)  +  (peteam-pchs(ib)+ccc2)*pzcom(ib) 

30  CONTINUE 
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Call  FINDMIN(cbar,ib,mini) 
cfmin  =  cbar(mini) 
oplihr  =  TL(mini) 
opuihr  =  TU(mini) 

TL2  =  oplihr 
TU2  =  opiithr 
F20  =  fusO 
F21  =  -fusO 
T2  =Tss 

RETURN 

END 


SUBROUTINE  SETBAND3(ccc3,TL3.TU3,F300.F301.F3l0,F31I.T31,T32) 
REAL  pO,  pi 

REAL  mhO,  mhl,  mslO,  msH, ms20,  tns21 
REAL  shO,  shl.  sslO,  ssll,  ss20,  ss2I 
REAL  cOO.clO.cll.cOl 
REAL  Tssl,  Tss2 

REAL  ccc3 

REAL  cfmin,  oplihr,  opuihr 
REALTU(lOOl),  TL(lOOl) 

REAL  pchs(lOOl),  pzcom(lOOl) 

REAL  cbar<  1001) 

p0  =  0.5 
pi  =0.5 


mhO  =  -1.0 
mhl  =  1.0 
mslO  =  -1.0 
msl  1  =  1.0 
ms20  =  -1.0 
ms21  =  1.0 

sigma  =  1.0 

cOO  =  0.0 
clO=  1.0 
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cl  1  =0.0 
cOl  =  1.0 

shO  =  sigma 
shl  =  sigma 
sslO  =  sigma 
ssl  1  =  sigma 
ss20  =  sigma 
ss21  =  sigma 

plamssl  =  (clO  -  c00)/(c01  -  cll) 
plamss2  =  (clO  -  c00)/(c01  -  cl  1) 
plamt  =  plamssl 

lO  =  pO/p  i 

Ths  =  (mhO  +  mhl)/2.0  +  (sigma**2/(mhl  -  mhO))*LOG(plami*lO) 
Tssl  =  (mslO  +  msll)/2.0 

•  +  (sigma**2/(msl  1  -  mslO))*LOG(plamssl*iO) 

Tss2  =  (ms20  +  ms2l)/2.0 

•  +  (sigma*  •2/(ms21  •  ms20))*LOG(piamss2*lO) 

bslO  =  (Tssl  -  mslOVsslO 
bsll  =  (Tssl  -  msll)/ssll 

CALL  QfunE(bslO,  QbslO) 

CALL  QfunE(bsll,Qbsll) 

bs20  =  (Tss2  -  ms20)/ss20 
bs21  =(Tss2-ms21)/ss21 

CALL  QfunE(bs20,  Qbs20) 

CALLQfunE(bs21,Qbs21) 

IDO  =  plamss  1  *  t0*(  1 .0-Qbs  1 0)*(  1 .0-Qbs20V((  1 .0-Qbs  1 1  )*(  1 .0-Qbs2 1  j) 
ID  1  =  plamss  1  •iO*(  1 .0-Qbs  1 0)*Qbs20/((  1 .0-Qbs  1  l)*Qbs2 1 ) 
no  =  plamss2*t0*Qbsl0*(1.0-Qbs20)/(Qbsll*(1.0-Qbs21)) 
fll  =  plamss2*i0*Qbsl0*Qbs20/(Qbsll*Qbs21) 

IDOhO  =  (IDO  -  mhOVshO 
fDlhO  =  (lDl  -  mhOVshO 
flOhO  =  (flO  -  mhO)/shO 
fllhO  =  (fll  -mhO)/shO 

CALL  QfunE((DOhO,  QlDOhO) 

CALL  Qfui\E(f01h0,  QlDlhO) 
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CALL  QfunE(flOhO,  QHOhO) 

CALL  QfunE(fllhO,  QfllhO) 

roohl  =  (fDO  -  mhD/shl 
fDlhl  =(fDl  -mhD/shl 
flOhl  =  (flO-  mhlVshl 
nihl  =(fll  -mhiyshl 

CALL  QfunE(fDOhl.  QfDOhl) 

CALL  QfunE(f01hl.  QrOlhl) 

CALL  QfunE(flOhl.  QnOhl) 

CALL  Qfui£(fllhl.Qnihl) 

peteam  =  (  Q(D0h0*(1.0-Qbsl0)*(1.0-Qbs20) 

•  +  QroihO*(l.O-Qbsl0)*Qbs2O 

•  +Qfl0h0*Qbsl0*(1.0-Qbs20) 

•  +  Qfl  lh0*Qbsl0*Qbs20 )  •  pO 

•  +(  (1.0-Qf00hl)*(1.0-Qbsll)*(1.0-Qbs2l) 

•  +(1.0-QfDlhl)*(1.0-Qbsll)*Qbs21 

•  +(1.0-Qfl0hl)*Qbsll*(1.0-Qbs21) 

•  +(1.0-QfllhirQbsll*Qbs21)*pl 

tine  =  0.02 
ib  =  0 

DO  30  ihr  =  Ths,  Ths  +  2.0,  tine 
ib  =  ib  +  1 
TU(ib)  =  ihr 

TL(ib)  =  Ths  -  (FLOAT(ib)  -  I.0)*tine 

tuhO  =  (TU(ib)  -  mhO)/shO 
tlhl  =  (TL(ib)  -  mhD/shl 

CALL  QfunE(tuhO.  QtuhO) 

CALLQfunE(tJhl,Qtlhl) 

pchs(ib)  =  QiuhO*pO  +  (1.0  -  QtJhl)*pl 

tlhO  =  (TL(ib)  -  mhOyshO 
tuhl  =  (TU(ib)  -  mhD/shl 

CALL  QfunE(UhO,  QdhO) 

CALL  QfunE(tuhl,QtuhD 

pzcom(ib)  =  (QiJhO  -  QtuhO)*pO  +  (Qtihl  -  QtuhD*pl 
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cbar(ib)  =  pehs(ib)  +  (peteam-pehs(ib)-H:cc3)*pzcom(ib) 


30  CONTINUE 

CALL  FINDMIN(cbar,ib.mini) 
cfmin  =  cbar(mini) 
oplihr  =  TL(mini) 
oputhr  =  TU(mini) 

TL3  =oplthr 
TU3  =  oputhr 
F300  =  fD0 
F301  =  roi 
F310  =  fl0 
F311  =  -fD0 
T31  =Tssl 
T32  =Tss2 

RETURN 

END 


SUBROUTINE  SETB AND23(ccc23.TL3 1  ,TL32,TU32.TU3 1 , 
•  F0E1.F00E01E10JF11,T231,T232) 


REAL  pO,  pi 

REAL  mhO,  mhl,  mslO,  msll,  ms20,  ms21 
REAL  shO,  shl,  sslO,  ssl  1,  ss20,  ss21 
REALcOO.clO.cll.cOl 
REAL  Tssl,  Tss2 

REAL  cccl,  ccc2 
REAL  cfmin 

REAL  optll,  optl2,  optu2,  opiul 
REAL  TUl(lOOl).  TU2(1001),  TLl(lOOI).  TL2(100I) 
REAL  pchs(lOOl),  pzcom  1(1001),  pzcom2(1001) 
REAL  cbar(lOOl) 

cccl  =  CCC23/2.0 
ccc2  =  ccc23 


p0  =  0.5 
pi  =0.5 
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mhO  =  -1.0 
mhl  =  1.0 
mslO  =  -1.0 
msll  =  1.0 
ms20  =  -1.0 
ms21=  1.0 

sigma  =  1.0 


cOO  =  0.0 
clO=  1.0 
cl  1=0.0 
cOl  =  1.0 

shO  =  sigma 
shl  =  sigma 
sslO  =  sigma 
ssll  =  sigma 
ss20  =  sigma 
ss21  =  sigma 

plamssl  =  (clO  -  c00)/(c01  -  cll) 
plamss2  a  (clO  -  c00)/(c01  -  cll) 
plamt  =s  plamssl 

lO  =pO/pl 

Ths  =  (mhO  +  mhl)/2.0  +  (sigma**2/(mhl  -  mhO))*LOG(plami*lO) 
Tssl  =  (mslO  +  msl  l)/2.0 

*  +  (sigiTia**2/(msll  -  mslO))*LOG(plamssl*iO) 

Tss2  =  (ms20  +  ms21)/2.0 

•  +  (sigma**2/(ms21  -  ms20))*LOG(plamss2*iO) 

bslO  =  (Tssl  -  mslO)/sslO 
bsll  =(Tssl  -  msll)/ssll 

CALL  QfunE(bslO,  QbslO) 

CALL  QfunE(bsIl,Qbsll) 

bs20  =  (Tss2  -  ms20)/ss20 
bs21  =  (Tss2  -  ms21)/ss21 

CALL  QfunE(bs20,  Qbs20) 

CALLQfunE(bs21,Qbs21) 

fsO  =  plamssl*i0*(1.0-Qbsl0)/(1.0-Qbsll) 
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fsl  =plamssl*tO*QbslO/Qbsll 

fsOhO  =  (fsO  -  mhOyshO 
fslhO  =  (fsl  -  mhOyshO 


CALL  QfunE(fsOhO,  QfsOhO) 
CALL  QfunE(fslhO.  QfslhO) 


fsOhl  =  (fsO  -  mhiyshl 
fslhl  =(fsl  -mhiyshl 

CALL  QfunE(fsOhl.  QfsOhl) 

CALL  QfunE(fslhl,  Qfslhl) 

peteaml  =  (Qfs0h0*(1.0-Qbsl0)  +  QfslhO*QbslO)  *  pO 
*  +  (( 1 .0-QfsOh  1  )•  ( 1 .0-Qbs  1 1 )  +  (1 .0-Qfs  Ih  1  )*Qbs  1 1 )  *  p  1 

roO  =  plamssl*t0*(1.0-Qbsl0)*(1.0-Qbs20y((1.0-Qbsll)*(1.0-Qbs21)) 
fDl  =  plamssl*t0*(1.0-Qbsl0)*Qbs20/((L0-Qbsn)*Qbs21) 
no  =  plamss2*t0*Qbsl0*(1.0-Qbs20y(Qbsl  I*(1.0-Qbs21)) 
fll  *  plamss2*t0*Qbsl0*Qbs20/(Qbsll*Qbs21) 


roOhO  =  (fOO  •  mhO)/shO 
fOlhO  =  ((01  -  mhOyshO 
flOhO  =  (flO  -  mhOyshO 
fllhO  =  (fll  -mhOyshO 


CALL  QfunE((OOhO.  QfOOhO) 
CALL  QfunE(roihO.  QfOlhO) 
CALL  QfunE(flOhO,  QnOhO) 
CALL  QfunE(f  1  IhO,  Qfl  IhO) 


(DOhl  =  (fOO  -  mhiyshl 
fOlhl  =  (f01  -  mhiyshl 
flOhl  =  (flO-  mhiyshl 
fllhl  =  (fll  -  mhiyshl 


CALL  QfunE((00hl.  QfOOhl) 
CALL  QfunE(f01hl.  QfDlhl) 
CALL  QfunE(flOhl,  QnOhl) 
CALL  QfunE(fllhl,Qnihl) 


peteam2  =  (  Qf00h0*(1.0-Qbsl0)*(1.0-Qbs20) 
"  +Q(01h0*(l. 0-Qbs  10)*Qbs20 

•  +Qfl0h0*Qbsl0*(1.0-Qbs20) 
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•  +  Qn  lh0*Qbsl0*Qbs20 )  *  pO 

•  +(  (l.0-QfD0hl)*(1.0-Qbsll)*(1.0-Qbs21) 

•  +(I.0-Qroihl)*(1.0-Qbsll)*Qbs21 

•  +  (1.0-Qfl0hl)*Qbsll*(1.0-Qbs21) 

•  +  (1.0-Qnihl)*Qbsll*Qbs21  )*pl 


tine  =  0.02 
ib  =  0 

EK)  30  ihr  =  Tbs,  Ths  +  2.0,  tine 
ib  =  ib  +  1 
TUl(ib)  =  thr 

TLl(ib)  =  Ths  -  (FLOAT(ib)  -  1.0)*tine 

TU2(ib)  *  (TUl(ib)+Ths)/2.0 
TL2(ib)  =  TLl(ib)+(Ths-TLl(ib))/2.0 

tulhO  =  (TUl(ib)  ■  mhOyshO 
tllhl  =  (TLl(ib)  -  mhiyshl 

CALL  QfunE(tulhO.  QtuIhO) 
CALLQfunE(aihl.QtUhl) 

pehs(ib)  =  QtulhO'pO  +  (1.0  -  Qtllhl)*pl 

UlhO  =  (TLl(ib)  -  mhOyshO 
tl2hO=»(TL2(ib)-mhOyshO 
tulhO  =  (TUl(ib)  -  tnhOyshO 
tu2h0  =  (TU2(ib)  -  mhOyshO 

CALL  QfunECUlhO,  QUlhO) 

CALL  QfunE(ll2hO.  QU2hO) 

CALL  QfunE(tulhO.  QtulhO) 

CALL  QfunE(tu2hO.  Qtu2hO) 

tllhl  =  (TLl(ib)- mhiyshl 
il2hl  =  (TL2(ib)- mhiyshl 
tulhl  =(TUl(ib)*  mhiyshl 
tu2hl  =(TU2(ib)- mhiyshl 

CALL  QfunE(aihl.Qillhl) 

CALL  QfunE(U2hl,Qil2hl) 

CALL  QfunE(tulhl,  Qlulhl) 

CALL  QfunE(tu2hl,Qlu2hl) 

pzcoml(ib)  =  (  QUlhO  -  Qll2hO  +  Qtu2hO  -  QtulhO )  *  pO 
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*  +  ( Qaihl  -  QiJ2hl  +  Qiu2hl  -  Qtulhl )  *  pi 
pzcom2(ib)  =  (Qtl2h0  -  Qlu2h0)*p0  +  (Qtl2hl  -  Qtu2hl)*pl 
cbar(ib)  =  pehs(ib) 

*  +(peteaml-pehs(ib)+cccl)*pzcoml(ib) 

*  +  (peteam2-pehs(ib)+ccc2)*pzcom2(ib) 

*  +  (pehs(ib)-peteam  1  -peifiam2-ccc  I  -ccc2) 

*  *pzcoml(ib)*pzcom2(ib) 

30  CONTINUE 

CALL  FINDMIN(cbar,ib,mini) 

cfmin  =  cbar(mini) 

opill  =  TLl(mini) 

optl2  =  TL2(mini) 

opiu2  =  TU2(mini) 

optul  =  TUI  (mini) 

TL31  =opai 
TL32  =  opU2 
TU32  =  optu2 
TU31  *  optul 
FO  =fs0 
FI  =  -fsO 
FOO  =f00 
FOl  =101 

Fio  =  no 
Fll  =-fOO 
T231  =  Tssl 
T232  =  Tss2 

RETURN 

END 


SUBROUTINE  QfunE(xx,  erfcx) 

REAL  X,  XX,  erfcx 

REAL  al,  a2.  a3,  a4,  a5,  p.  pi,  t 

pi  =  3.141592654 
X  =  ABS(xx) 


al  =  0.319381530 
a2  =  -0.356563782 
a3  =  1.781477937 
a4  =  -1.821255978 
a5=  1.330274429 

p  =  0.2316419 
t=  1.0/(1.0-t-p*x) 


si  =  al*t  +  a2*i**2  -►  a3*t**3  +  a4*t**4  a5*t**5 
s2  =  sl*EXP(-(x‘*2)/2.0) 

IF  (XX  .GE.  0.0)  THEN 
erfcx  *  s2/SQRT(2.0*pi) 

ELSE  IF  (XX  .LT.  0.0)  THEN 
erfcx  =  1.0  -  s2/SQRT(2.0*pi) 

END  IF 

RETURN 

END 


0%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

SUBROUTINE  FINDMIN(array,isize,minindex) 

REAL  array(isize) 

INTEGER  minindex 

int  =  1 

11  CONTINUE 

DO  10  i  =  int+1,  isize 
IF  (arTay(int)  .GT.  aiTay(i))  THEN 
minindex  =  i 
int  =  minindex 
GOTO  11 
END  IF 

10  CONTINUE 

RETURN 

END 
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